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Motor  calculus  is  a  mathematical  system  that  is  quite 
analogous  to  vector  calculus  and  is  particularly  well-suited 
for  rigid  body  dynamics.  This  work  presents  a  new  approach 
to  motor  calculus  and  its  application  to  dynamics.  As 
opposed  to  the  common  bivector  definition,  a  motor  is  here 
defined  to  be  a  special  kind  of  vector  field.  It  is 
demonstrated  that  this  novel  motor  definition  is  both  useful 
for  investigating  the  mathematical  structure  of  motor 
calculus  and  for  applying  motor  calculus  to  rigid  body 
dynamics.  New  results  include  simple  derivations  and 
representations  of  the  inertia  dyadic  for  a  rigid  body  and 
of  the  equations  of  motion  for  a  serially  connected  open- 
loop  kinematic  chain. 


CHAPTER  I 
INTRODUCTION 

Motor  calculus  is  an  elegant  mathematical  system  that 
is  particularly  well-suited  to  rigid  body  kinematics  and 
dynamics.  Briefly,  a  motor  is  a  geometric  entity  that  can 
be  identified  with  the  following  physical  entities:  the 
velocity  distribution  of  a  rigid  body,  the  moment 
distribution  of  a  force  system,  and  the  momentum 
distribution  of  a  mass  system.  These  vector  fields  are 
respectively  called  the  velocity  motor,  the  wrench  and  the 
momentum  motor.  Rigid  body  dynamics  can  be  formulated  as 
the  study  of  the  algebraic  and  differential  relationships 
among  these  three  kinds  of  vector  fields;  these 
relationships  are  studied  by  means  of  motor  calculus. 

The  basic  ideas  of  motor  calculus  received  considerable 
attention  in  the  late  19th  century.  Perhaps  the  single  most 
important  work  of  that  time  was  a  treatise  by  R.  S.  Ball 
entitled  The  Theory  of  Screws  (1900).  A  screw  is  defined  as 
a  line  together  with  a  real  number  called  the  pitch.  By 
adjoining  another  real  number,  called  the  magnitude,  to  the 
screw,  we  obtain  an  entity  that  can  be  identified  with  a 
motor.  It  was  Clifford  (1873),  in  fact,  who  defined  the 
motor  as  a  "magnitude  associated  with  a  screw." 


Ball  considered,  as  did  many  other  authors  at  that  time 
(perhaps  most  notably  Plucker,  1866),  the  dynamics  of  a 
rigid  body  from  the  standpoint  of  line  geometry;  the  screw 
is  considered  the  fundamental  entity.  The  magnitude 
associated  with  a  screw  was  incorporated  through  the 
coordinates  used  to  represent  the  screw — the  screw 
coordinates  as  Ball  called  them.  Ball  derived  many  of  the 
dynamical  relations  of  an  initially  quiescent  body  in  terms 
of  screw  coordinates  (he,  in  effect,  did  not  consider 
inertial  forces  and  moments  due  to  velocity). 

In  motor  calculus  the  motor  as  opposed  to  the  screw  is 
considered  to  be  the  fundamental  entity.  By  introducing 
various  algebraic  and  differential  operations  for  motors, 
von  Mises  (1924a  and  1924b)  developed  motor  calculus  in  such 
a  way  that  it  is  quite  analogous  to  vector  calculus.  He 
also  illustrated  the  ease  with  which  motor  calculus  can  be 
applied  to  the  dynamics  of  a  rigid  body  in  general  (he 
included  the  inertial  forces  and  moments  due  to  velocity). 

For  some  forty  years  after  von  Mises 's  works,  screw 
theory  and  motor  calculus  received  little  attention.  The 
revival  of  these  subjects,  much  of  which  has  been  in  the 
area  of  instantaneous  kinematics  of  rigid  bodies,  was 
initiated  by  Phillips  and  Hunt  (1964),  and  there  is 
continuing  interest  at  the  present  time  (for  a  good  review 
of  recent  contributions,  see  Hunt,  1978;  Bottema  and  Roth, 
1979;  and  Phillips,  1984).   Significant  recent  developments 


in  motor  calculus  and  its  application  to  rigid  body  dynamics 
are  due  to  Dimentberg  (1965),  Yang  (1967,  1969,  1971,  1974), 
Woo  and  Freudenstein  (1971),  Pennock  and  Yang  (1982),  and 
Featherstone  (1983,  1984). 

It  was  my  original  intention  to  consider  the  use  of 
motor  calculus  in  developing  efficient  algorithms  for 
dynamic  analysis  of  serially  connected  open-loop  kinematic 
chains  ( it  should  be  pointed  out  that  this  has  been 
investigated  most  recently  by  Featherstone,  1984).  It  turns 
out,  however,  that  the  emphasis  of  this  work  is  more 
theoretical  than  originally  intended.  In  the  course  of  the 
research  it  became  apparent  that  there  were  many  areas  of 
motor  calculus  proper  and  its  application  to  rigid  body 
dynamics  that  could  be  treated  more  rigorously  than  had  been 
done  previously. 

This  dissertation  is  essentially  a  rigorous  development 
of  motor  calculus  and  its  application  to  rigid  body 
dynamics.  An  important  new  concept  introduced  is  the 
definition  of  a  motor  as  a  special  kind  of  vector  field. 
This  is  in  contrast  to  the  common  bivector  definition  (here, 
we  use  the  bivector  as  a  particular  representation  for  a 
motor).  The  motivation  for  the  vector  field  definition  is 
that  it  is  more  tractable  as  regards  analytical 
considerations  of  motor  calculus  and  as  regards  the 
application  of  motor  calculus  to  dynamics.  Aside  from 
investigating    motor    calculus    proper    and    describing 


elementary  dynamics  via  motor  calculus,  we  also  consider  the 
dynamics  of  serially  connected  open-loop  kinematic  chains 
using  some  of  the  new  concepts  developed. 

We  now  briefly  describe  the  contents  each  of  the 
chapters . 

In  Chapter  II  we  define  the  motor  as  a  vector  field 
with  a  special  property,  and  we  then  proceed  to  develop  the 
fundamentals  of  motor  calculus.  With  the  new  motor 
definition,  addition  of  motors  is  defined  simply  as 
pointwise  addition  of  the  corresponding  vector  fields  and  we 
show  that  the  derivative  of  a  motor  is  given  by  pointwise 
differentiation  of  the  vectors  in  the  corresponding  vector 
field.  It  is  believed  that  some  of  the  results  obtained  for 
motor  dyadics  and  motor  derivatives  and  integrals  are  new. 
The  dyad  expansion  that  is  given  for  a  motor  dyadic  will  be 
seen  to  have  novel  and  important  applications  in  the 
formulation  of  the  dynamic  equations  for  a  rigid  body,  and 
it  is  shown  that  by  considering  the  set  of  all  motors  as  a 
normed  space  enables  limiting  operations  for  motor  functions 
to  be  treated  in  a  more  general  and  rigorous  manner. 

In  Chapter  III  we  consider  various  geometrical 
properties  of  motors  and  motor  operations.  We,  in  effect, 
consider  some  of  the  connections  between  motor  calculus  and 
screw  theory.  These  geometrical  aspects  of  motor  calculus 
are  often  useful  when  motor  calculus  is  applied  to  rigid 
body  kinematics  and  dynamics.    Aside  from  using  the  new 


motor  definition  to  describe  geometric  properties,  a  novel 
derivation  of  the  so-called  circle  representation  of  a  two- 
system  is  given,  which  is  believed  to  be  simpler  than  any 
other  found  in  the  literature. 

In  Chapter  IV  we  consider  elementary  dynamics  using  the 
motor  calculus  developed  in  Chapter  II.  Here,  the  new  motor 
definition  is  particularly  useful  for  casting  the  laws  of 
motion  in  terms  of  motors. 

In  Chapter  V  we  consider  elementary  kinetostatics  using 
material  from  both  Chapters  II  and  III.  We  consider  in 
particular  the  static  and  velocity/acceleration  analysis  of 
an  open-loop  chain.  The  static  analysis  for  an  open-loop 
chain  constitutes  a  generalization  of  that  found  in  Lipkin 
and  Duffy  (1982)  and  the  equations  given  here  for 
velocity/acceleration  analysis  of  an  open-loop  chain  are  the 
same  as  those  given  in  Featherstone  (1984). 

In  Chapter  VI  we  tie  together  many  of  the  results  of 
Chapters  II-V  for  application  to  rigid  body  dynamics.  The 
new  definition  of  a  motor  enables  a  relatively  simple 
construction  of  the  inertia  dyadic  for  a  rigid  body. 
Further,  the  dyad  expansion  representation  of  the  inertia 
dyadic  given  here  is  novel  and  for  the  first  time  it  relates 
directly  to  the  classical  principal  screws  of  inertia  due  to 
Ball  (1900).  Finally,  the  equations  of  motion  for  an  open- 
loop  chain  are  derived.  As  far  as  the  author  is  aware,  both 
the  derivation  and  the  form  of  these  equations  are  the 
simplest  to  appear  in  the  literature  thus  far. 


CHAPTER  II 
FUNDAMENTALS  OF  MOTOR  CALCULUS 

In  this  chapter  motor  calculus  is  developed  using  the 

novel  vector  field  definition  of  the  motor.    We  assume 

familiarity  with  elementary  vector  calculus,  linear  algebra, 

and  calculus. 

2.1   The  Motor 


We  denote  vectors  and  vector  quantities  by  bold  lower 
case  letters  (e.g.,  v) ,  and,  as  is  customary  in  engineering, 
we  shall  usually  confound  these  two  terms.  The  association 
or  lack  thereof  of  a  vector  with  a  point  or  line  determines 
whether  the  vector  is  considered  free  or  bound.  Whether  a 
vector  is  being  considered  a  free,  point,  or  line  vector 
will  often  not  be  made  explicit  (this  is  also  common  in 
engineering).  In  fact,  the  same  vector  may  have  different 
associations  in  different  contexts.  The  reason  for  this 
convention  is  that  the  association  (or  lack  thereof)  is 
often  clear  without  making  formal  algebraic  distinctions. 
For  example,  suppose  that  vectors  u  and  v  are  associated 
with  two  distinct  points.  Then,  we  assume  that  such  an 
expression  as  u  =  v  implies  that  u  and  v  have  the  same 
magnitude  and  direction.  Or,  the  sum  u  +  v  is  obtained  by 
assuming  u  and  v  are  associated  with  the  same  point  (i.e.. 


assuming  that  u  and  v  are  translated  in  a  parallel  sense  to 
tne  same  point  so  that  they  can  be  added),  and  then  this  sum 
may,  depending  on  the  context,  be  a  free,  point,  or  line 
vector. 

A  vector  field  is  a  set  of  vectors  that  is  indexed  by 
all  points  in  space.  A  vector  at  some  point  P  in  a  vector 
field  v  is  denoted  Vp .  Here,  v  can  be  interpreted  as  an 
arbitrary  vector  in  the  vector  field  to  which  it  belongs. 
We  now  define  the  motor  as  a  special  kind  of  vector  field 
and  give  some  associated  definitions  and  properties. 

DEFINITION.  A  motor  is  a  vector  field  m  for  which 
there  is  a  vector  in  such  that  for  every  pair  of  points  P,  Q, 
we  have 

mQ-mp  =  inxPO,  (2.1) 

where  PQ  is  a  vector  directed  from  P  to  Q. 

The  vector  m  is  called  the  principal  vector*  of  the 
motor,  and,  as  regards  the  above  definition,  this  vector  is 
considered  free.  In  mechanics,  however,  it  is  often 
convenient  to  consider  the  principal  vector  bound. 

THEOREM  2.1.  The  principal  vector  of  a  motor  is 
unique . 


*    This  term  is  borrowed  from  Dimentberg  (1965)  and  its 
meaning  here  is  equivalent  to  his  definition. 


PROOF.  Suppose  in  and  in'  are  principal  vectors  of  the 
motor  m  and  let  P,  Q  be  chosen  as  distinct  points  such  that 
PQ  •  in  =  PQ  •  in'  =0  (i.e.,  PQ  is  perpendicular  to  m  and 
m' ) .   From  eqn.  (2.1),  we  must  have 

m  X  PQ  =  m'  X  PQ    , 

or,  on  rearranging, 

(in  -  m')  X  PQ  =  0 

Cross  multiplying  both  sides  with  PQ  and  using  the  identity 
a  X  (b  X  c)  =  (a  •  c)b  -  (a  •  b)c,  we  obtain 

(PQ  .  PQ)(m  -  m'  )  -  [PQ  •  (m  -  in')]PQ  =  0    .   (2.2) 

Since  PQ  /^  0  ( P  and  Q  are  distinct  by  assumption)  and 

PQ    •    (fii  -   m' )    =    PQ    •    m   -    PQ    •   m' 

=    0    +    0 

eqn.    (2.2)    becomes 

iff   -    nf'    =    0 


Thus  ra  =  m' .  Q.E.D. 

We  denote  motors  by  bold  capital  letters  (e.g.,  M,  N) 
and  the  set  of  all  motors  by  _M.  The  notation  M:m  or  m:M 
indicates  the  vector  field  m  corresponding  to  M.  When  we 
write  upper  case  letters  for  motors  (M,  N,  L,  etc.),  we 
shall  always  assume  that  vectors  in  their  corresponding 
vector  fields  are  denoted  by  corresponding  lower  case 
letters  (m,  n,  1,  etc.).  Thus,  if  the  vector  field  v  is  a 
motor,  and  we  wish  to  denote  this  motor  by  M  (i.e.,  M:v)  , 
then  we  have  m  =  v  (where  we  assume  that  m  and  v  are 
associated  with  the  same  point). 

Two  motors  M,  N  are  said  to  be  equal,  M  =  N,  if  and 
only  if  m  =  n  (i.e.,  the  vector  fields  are  equal  pointwise). 
It  follows  from  Theorem  2.1  that  m  =  n. 

Finally,  the  zero  motor  M  =  0  is  defined  as  the  zero 
vector  field  m  =  0.   It  follows  from  eqn.  (2.1)  that  in  =  0 . 


2. 2   The  Bivector  Representation 


DEFINITION.  A  bivector  corresponding  to  a  motor  M  is 
the  ordered  pair  of  vectors  (m,  rap)  where  P  can  be  any 
point. 

THEOREM  2.2.  A  motor  is  characterized  by  any  one  of 
its  bivectors  and  any  bivector  characterizes  a  unique  motor. 

PROOF.  This  follows  immediately  trom  eqn.  (2.1)  and 
Theorem  2.1.   Q.E.D. 
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By  virtue  ot  this  theorem,  any  motor  can  be  represented 
by  a  necessarily  unique  bivector  at  a  given  point.  We 
denote  the  set  of  all  bivectors  at  a  point  P  by  ^p. 
Clearly,  M_  and  _Bp  have  the  same  cardinal  number.*  We 
indicate  the  bivector  of  a  motor  M  at  a  point  P  by  (M)p; 
thus,  (M) p  =  (m,  mp) . 

Two  bivectors  (M)p,  (N)p  are  said  to  be  equal,  (M)p  = 
(N)p,  if  and  only  if  m  =  n  and  mp  =  np .  It  follows  that 
(M)p  =  (N)p  if  and  only  if  M  =  N. 

Finally,  the  zero  bivector  is  defined  as  (0,0).  It 
follows  that  (0)p  =  (0,0). 

2 . 3  The  Velocity  Motor  and  the  Wrench 
We  introduce  here  two  of  the  three  kinds  of  motors  with 
which  we  shall  be  concerned  in  rigid  body  mechanics.  As 
mentioned  in  the  introduction,  applying  motor  calculus  to 
rigid  body  dynamics  comprises  the  study  of  the  algebraic 
and  differential  relationships  among  wrenches,  velocity 
motors  and  momentum  motors.  We  introduce  the  velocity 
motor  and  wrench  here  so  that  we  can  discuss  some  of  their 
mathematical  properties  in  subsequent  sections  of  this 
chapter  and  in  the  next  chapter.  We  reserve  the 
introduction  of  the  momentum  motor  for  Chapter  IV.    (It 


That  is,  there  exists  a  one-one  correspondence  between  M 


and  tir,. 
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turns  out  that  the  definition  of  the  momentum  motor  is  quite 
similar  to  that  of  the  wrench.) 

DEFINITION.  The  velocity  distribution  of  one  frame 
relative  to  another  is  called  a  velocity  motor. 

Let  V  be  the  velocity  distribution  of  frame  S'  relative 
to  frame  Z.  That  v  is  in  fact  a  motor  is  verified  by 
differentiating  (with  respect  to  time)  a  vector  joining  two 
fixed  points  P,Q  in  E' .   On  the  one  hand  we  have 


1^  PQ  =  CO  X  PQ 


where  o)  is  the  so-called  angular  velocity  of  Z'  relative  to 
T.,    and  on  the  other  hand 


it  PQ  =  ^Q  -  ^P 


Thus, 


Vq  -  Vp  =  0)  X  PQ 


The  velocity  motor  V:v  has  principal  vector  Q'  =  u. 

DEFINITION.   The  moment  distribution  of  a  force  system 
is  called  a  wrench.* 


*   Ball  (1900)  also  used  the  term  wrench  but  not  in  the  same 
sense  as  here.   We  consider  Ball's  definition  in  Chapter  IV. 
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Let  m  be  the  moment  distribution  of  a  force  system, 
composed  of  forces  f  ^^ ,  ±2,  ...,  f^  (these  force  vectors  can 
be  considered  point  vectors).  In  order  to  show  that  ni  is  a 
motor,  let  P,Q  be  arbitrary  points  and  let  r^,  r^  be  vectors 
directed  from  P  and  Q  to  the  point  of  application  of  f  •  . 
Then,  we  have 


">Q  -  "p  =  I  rj  X  f^  -  J  r^  X  f . 


=  Krj  -  4,  ,  .. 


=  I  QP  X  f . 


=  (I  f^)  X  PQ 


Thus,  the  wrench  M:m  has  principal  vector  m  =  If.   . 

Because  of  the  nature  of  the  cross  product,  the  forces 
can  be  slid  along  their  lines  of  action  without  altering  the 
moment  distribution.  Thus,  if  desired,  the  forces  of  a 
force  system  can  be  considered  line  vectors  as  regards  the 
wrench  of  this  force  system. 

We,  of  course,  can  be  more  general  and  say  that  the 
moment  distribution  of  any  set  of  point  (or  line)  vectors  is 
a  motor.  We  shall  use  this  fact  when  we  consider  the 
momentum  motor  in  Chapter  IV. 
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2 . 4   Algebraic  Operations 
For  every  pair  of  motors  M,  N  and  real  scalars  a,  6 
(or,  more  concisely,   for  all  M,  N  e   M_  and  a,  B  e  ^)  ,   we 
define  scalar  multiplication  and  add  it  ion/subtract  ion  by 

otM  ±  6N  :  am  ±  6n    ; 

the    cross   product    by 

MxN:mxn   +   mxn         ; 

and  the  reciprocal  product  by 


MoN  =  iii«n+ia^  .n 


where  P  is  any  point.    In  addition,  we  define  -M:-m  (or, 
equivalently ,  -M  =  (-l)M). 

In  order  to  verify  that  the  vector  fields  am  ±  5n  and 
m  X  n  +  m  X  n  are  in  fact  motors,  we  must  show  that  they 
have  principal  vectors  (such  that  eqn .  (2.1)  is  satisfied). 
For  every  pair  of  points  P,0  we  have 


( omp,  ±  BOp.)  -  (amp  ±  BOp)  =  a(mn  -  mp)  ±  6(0^  -  Op) 


Q    "P' 


•Q    "P' 


=  am  X  PQ  ±  Bn  X  PQ 


=  (  am  ±  Bn)  ^    PQ 
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and 


(m    X   n      +   m      X    n)    - 


(m   X   Dp   +  nip    X   n)    =  m   x    (ng  -    Hp)    +( nig  -  nip)     x  n 


=   inx    (nxPQ)    +    (mxPQ)     xm 

=    (m    X   n)     X    PQ         , 

which  verify  that   am  ±   6m  and  in  x  fi  are  the  principal 
vectors  of  aM  ±  BN  and  M  x  n  respectively. 

The  reciprocal  product  M  o  n  is  independent  of  the 
choice  for  P  since  for  every  point  Q,  it  can  be  shown  using 
eqn.  (2.1)  that 

m*nQ  +  mQ*m  =  in*np  +  mp*n 

We  note  that  since  M_  is  closed  under  scalar 
multiplication  and  addition  (i.e.,  aM  +  SN  e  _M_  for  all  M, 
N  e  M_  and  a,  B  e  ^)  ,  _M  is  a  vector  space.* 

We  now  define  analogous  operations  to  those  above  on 
the   set   of   bivectors  _Bp  such   that  an   ismorphism**   is 


*   That  all  the  axioms  defining  a  vector  space  are  satisfied 
is  easily  verified. 

**   An   isomorphism   is   a   one-one   correspondence   that 
preserves  significant  relations  between  two  sets. 
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established  between  _Bp  and  _M_.    For  all  M,  N  e  _M^  and  a, 
3  e  _R ,  the  operations  are  defined  by 

ct{m,  fflp)  ±  3(n,  Dp)  =  (  am  ±  3n,  amp  +  Bnp)    , 

(m,  mp)  X  (n,  np)  =  (m  x  n,  in  x  np  +  mp  x  FT)    , 

(m,  mp)  o  (n,  Op)  =  m  •  np  +  mp  •  n 

Then,  we  have  an  isomorphism  for  which 

(aM  ±  BN)p  =  a{M)p  ±  8(N)d 


(M)p  X  (N)p  =  (M  X  N)p 


(M)p  o  (N)p  =  M  o  N 


Note  that  ^  is  a  vector  space  and  that  ^'^  P  *—  *  -^  ^^  ^" 
onto,  invertible  linear  transformation. 

Finally,  we  define   (•)p   ^s  the  inverse  of  (•)„. 

It  is  interesting  to  note  that  with  the  above 
definition  of  motor  addition,  we  can  express  a  wrench  as  a 
sum  of  motors  each  of  which  is  the  moment  distribution  of  a 
single  force.  This  then  would  provide  an  alternative  method 
for  verifying  that  the  moment  distribution  of  a  force  system 
is  a  motor.   We  need  only  prove  that  the  moment  distribution 


1€ 


of  a  single  force  is  a  motor.   The  sum  then  is  a  motor  since 
M_  is  closed  under  addition. 

2 . 5   Some  Linear  Algebra  Concepts 
2.5.1   Coordinates  and  Subspaces 

THEOREM  2.3.   dim  _M  =  dim  ^  =  6.* 

PROOF.  Since  there  exists  an  invertible  linear 
transformation  of  _M_  onto  ^,  namely,  (Op/  we  have  that 
dim  _M   =   dim  B^.  Now,   let  {e  ,  e  ,  e  }  be   a   linearly 

independent   set   of   vectors.     Then,   the   set  {(e  ,  0), 
(62,  0),  (03,  0),  (0,  e^),  (0,  62),  (0,  e^)}  is  clearly  a 
basis  for  ^.   Thus  dim  _Bp  =  6.   Q.E.D. 

An  isomorphism  (with  respect  to  addition  and  scalar 

multiplication)  between  _M  and  B^    can  be  established  once  a 

basis  is  selected  for  M.   Let  {M,  ,  M_  ,  ...,  M^ }  be  a  basis 

for   M_.     Then   the   scalars   a^^ ,   02^   ...,   ag   such   that 

6 
M  =  I      a . M.  are  called  the  coordinates  of  M  with  respect  to 

i  =  l 
this  basis.    We  use  the  common  notation  [M]   to  denote  a 

column  vector  of  coordinates  (i.e.,  a  coordinate  vector): 


dim  M_  and  dim  _Bp  refer  to  the  dimensions  of  these  vector 

spaces;  the  dimension  of  a  vector  space  is  the  number  of 
elements  in  a  basis. 
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[M]  = 


'°1  ' 

°2 

«3 

«4 

°5 

.  "6  . 

By  assuming  the  standard  definitions  for  addition  and  scalar 

multiplication  of  matrices  (here,  [M]  is  a  6  x  1  matrix),  we 

establish   an   isomorphism   between   _M   and  R^ .  Since 

6 
(Mg)p  =  I       ot.(M.)p   and    {(M^)p,  (M)  ,  ...,  (M)  }   is  a 

i  =  l 
basis  for  _Bp ,  we  have  that  a^,     a.2t      •'•,     cig  ace  also  the 

coordinates  of  (M)p  with  respect  to  this  basis  of  bivectors. 

Let  E-^,      E2  f   ...f   Eg  be  defined  such  that   (  E-i  )  p  = 

(e^,  0),  (E2)p  =  (62,  0),  (E3)p  =  (63,  0),  ( E4 )  p  =  (0,  e^^), 

^^5^P  ~  ^^'  ®2^'  ^^6^P  ~  ^^'  ®3^  where  {e  ,  e  ,  e  }  is  an 

orthonormal  set  of  vectors.    Then,  {E. ,  E^ ,  ...,  E^ }  is  a 

12         6 

basis,  which  we  call  a  standard  basis  at  point  P.  Such  a 
basis,  it  turns  out,  is  often  convenient  to  use  and  is 
commonly  used  in  the  literature.*  When  we  refer  a  motor  M 
to  a  standard  basis  we  shall  often  simply  write 


*  Featherstone  (1984)  defined  the  term  standard  basis 
equivalently.  We  shall  show  in  Chapter  III  that  there  is  a 
rectangular  set  of  coordinate  axes  associated  with  every 
standard  basis. 
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[M]  = 


m 


m. 


instead  of  writing  the  six  coordinates. 

We  denote  the  subspace  spanned  by  any  set  _S_  of  motors 
by  <^> .  If  _N  is  an  n-dimensional  subspace  of  jA_  and 
{M,  ,  M„,  ...,  M  }  is  a  basis  for  N_,  then  we  have  that  <M-,  , 
M2  /  ...f  N  >  =  J^  where  for  simplicity  we  have  dropped  the 
set  brackets.  We  define  the  coordinates  of  a  motor  M  e  J^ 
with  respect  to  the  basis  for  _N_  in  a  like  manner  as  above; 
thus, 


t»*n^  = 


"1 
«2 


n 
where   M  =   7   a.  M.  . 
i  =  l  '    ' 


On   selecting   a  basis   for  _N,   an 


isomorphism  is  established  between  _N  and  ^ 


2.5.2   On  the  Reciprocal  Product 

The  reciprocal  product  is  a  special  kind  of  bilinear 
scalar  product  on  _M_.  It  is  particularly  useful  in  rigid 
body  dynamics  (  in  ways  analogous  to  the  usefulness  of  the 
dot  product  in  particle  dynamics).  Here,  we  give  some 
definitions  and  properties  concerning  this  product. 
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DEFINITIONS.  Two  motors,  M,  N  are  said  to  be 
reciprocal  to  each  other  (or,  simply,  reciprocal)  if  M  o  n  = 
0.   A  motor  M  is  sel f-reciprocal  if  M  o  m  =  0. 

Unlike  the  dot  product  for  vectors,  the  reciprocal 
product  is  not  positive  definite;  that  is,  there  exist 
motors  that  are  self-reciprocal.  For  example,  if  (M)p  = 
(m,  0)  or  (M)p  =  (0,  m)  ,  then  M  o  M  =  0.  The  reciprocal 
product,  however,  is  non-degenerate;  that  is,  when  the 
reciprocal  product  is  expressed  in  matrix  form, 

M  o  N  =  [M]'^Q[N1 

the  6x6  matrix  Q  has  full  rank  (or  equivalently ,  is  non- 
singular)  for  all  bases  that  induce  this  matrix  form.  Now, 
let  Q^-  denote  the  i,j^^  entry  of  Q.   Then, 


0-   =  M.  o  M. 


where  {M  ,  M  ,  ...,  M  }  is  a  basis  to  which  the  coordinates 

of  [M] ,  [N]  refer.  Referring  to  a  standard  basis,  we  have 

'o   0  0   1   0   0 

0   0  0   0   10 

0   0  0   0   0   1 

10  0   0   0   0 

0   10   0   0   0 

0   0  10   0   0 
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for  which  |q|  =  -1  (thus  verifying  that  Q  has  full  rank).  A 
diagonal  form  for  Q  is  obtained  by  referring  to  a 
co-reciprocal  basis;  that  is,  a  basis  for  which 


M^  o  M.  =  0      for      i  *  J 


Let  {e^,  e^ ,  e  }  be  an  orthonormal  set.   Then,  the  set  {M  , 
M2 /  .../  M  }  such  that 


(M^)p  =  (e^,  e^),  (H^)^   =    (e^,  -e^). 


(M3)p  =  (e^,  e^),  ( M^ ) p  =  ( e^ ,  -e^). 


(M5)p  =  (03,  63),  (M^)p  =  (e3,  -e3) 


is    a   particular    co-reciprocal    basis    for    which 


Q  = 


2  0  0  0  0  0 
0-20000 
0  0  2  0  0  0 
0  0  0-200 
0  0  0  0  2  0 
0       0       0       0       0-2 
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In  this  diagonal  form,  we  see  that  the  signature*  of  Q  is 
zero  and  that,  again,  Q  has  full  rank. 

Ball  (1900)  worked  almost  exclusively  with  co- 
reciprocal  bases.  (In  fact,  the  terms  reciprocal  and 
co-reciprocal  are  due  to  Ball.)  The  screw  coordinates  that 
Ball  used  are  equivalent  to  the  coordinates  of  a  motor  that 
is  referred  to  a  co-reciprocal  basis.  Most  of  the  equations 
in  Ball's  work  are  expressed  in  terms  of  these  coordinates. 
It  must  be  emphasized  that  it  is  considered  desirable  in 
this  dissertation  to  obtain  equations  and  expressions  in  the 
more  general  coordinate-free  form.  A  special  kind  of  motor 
dyadic  defined  by  von  Mises  (1924a)  (which  is  covered  in  the 
next  subsection)  enables  an  elegant  means  of  expressing  many 
of  Ball's  results  concerning  the  dynamics  of  a  rigid  body  in 
coordinate-free  form.  The  importance  of  the  reciprocal 
product  and  co-reciprocal  bases  is  immediately  manifest  in 
this  dyadic  (as  we  shall  see). 

As  regards  the  co-reciprocal  basis  given  above,  which 
was  used  to  obtain  the  diagonal  form  for  Q,  we  note  that 
each  motor  in  this  basis  is  not  sel  f -reciprocal .  We  can 
make  a  general  statement  concerning  this. 


*  The  signature  of  a  bilinear  form  is  defined  as  the  number 
of  positive  eigenvalues  minus  the  number  of  negative 
eigenvalues  of  its  representing  matrix  (for  our  case  Q)  . 
The  signature,  like  the  rank,  is  independent  of  the  choice 
of  basis. 
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THEOREM  2.4.  Let  the  set  of  non-zero  motors  {M,  ,  M2  / 
...,  Mg}  be  co-reciprocal.  Then,  this  set  is  linearly 
independent  if  and  only  if  each  motor  in  the  set  is  not 
self -reciprocal . 

PROOF.  We  must  show  that  a-,  =  oo  =  . .  .  =  otc  =  0  is  the 
only  solution  for  ot-j^M-j^  +  a2M2  +  ...  +  <^6^6  ~  ^*  f^ow ,  for 
some  j , 


J      11     22         66      3 


or. 


a.M.  o  M.  =  0 
2    J     J 


If  M^  is  not  self-reciprocal,  then  it  follows  that  a^  =  0. 

Conversely,  suppose  that  M-  is  self-reciprocal.  Assume 
that  the  set  (M^  ,  M_ ,  ...,  M^ }  is  linearly  independent. 
Now,  if  we  choose  this  basis  to  obtain  the  matrix  expression 
for  the  reciprocal  product,  that  is  [M]'^0[N],  we  find  that 
0,  which  is  a  diagonal  matrix,  has  a  zero  diagonal  entry 
corresponding  to  M-  o  m^ .  This  implies  that  0  is  less  than 
full  rank,  which  contradicts  the  assumption  of  linear 
independence.   Q.E.D. 

(Note  that  this  theorem  can  be  generalized  to  any  non- 
degenerate  scalar  product  that  is  not  positive  definite.) 
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DEFINITIONS.  Let  _N  c  M_  be  a  subspace .  Then,  the 
reciprocal  subspace  ^  is  the  set  of  all  motors  reciprocal 
to  every  motor  in  _N_.  The  two  subspaces  U_,  _N  are  said  to  be 
reciprocal  to  each  other  (or,  simply  reciprocal).  The 
subspace  _N_  is  said  to  be  sel f -reciprocal  if  ^  <=  ^^. 

That  _N_  is  a  subspace  follows  from  the  bilinearity  of 
the  reciprocal  product. 

Since  the  reciprocal  product  is  a  special  kind  of  non- 
degenerate  scalar  product,  theorems  concerning  these 
products  in  general  apply  to  the  reciprocal  product 
inparticular .  We  state  two  such  theorems  whose  proofs  can 
be  found  in  (Lang,  1966).* 

THEOREM  2.5.  Let  _N_^  be  reciprocal  to  _N_.  Then, 
dim  _N^  +  dim  _N  =  6. 

THEOREM  2.6.  Every  subspace  of  _M_  has  a  co-reciprocal 
basis . 

In  Lang,  these  theorems  result  from  a  rather  general 
treatment  of  non-degenerate  scalar  products.  A  less 
abstract  approach  can  be  found  in  (Sugimoto  and  Duffy, 
1982),  where  they  exploit  a  one-one  correspondence  between 
so-called  reciprocal  and  orthogonal  screw  systems  (these 
screw  systems  are  what  we  would  call  subspaces). 


see  Lang's  Theorem  5  on  p.  129  and  Theorem  9  on  p.  135. 
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It  is  important  to  note  that  a  subspace  J>]_  and  its 
reciprocal  subspace  _N^  need  not  be  complements;*  this  is 
because  of  the  existence  of  sel t-reciprocal  motors. 
Consider,  for  example,  the  subspace  _N  =  <M>  where  M  is  self- 
reciprocal.  It  then  follows  that  _N  is  self-reciprocal;  that 
is,  J;J_  c  _N^.  Thus  _N  and  ti^  are  not  complements  since  _N  n  N^ 
=  ^  ^  {0}.  Finally,  note  that  a  self-reciprocal  subspace 
can  have  dimension  no  greater  than  three.  For,  if  N  is 
self-reciprocal,  then  we  must  have  dim  N_  <  dim  _N^  (as 
_N  c  _N_^)  ,  and  Theorem  2.5  would  be  contradicted  it  dim  J;J_  >  3. 

2.5.3   Dyadics,    Induced   Inner   Products   and   Norms   of 
Symmetric  Dyadics 

We  represent  linear  transformations  on  _M  by  means  of 
dyadics.  Motor  dyadics  were  originally  employed  by 
von  Mises  {1924a  and  1924b)  and,  as  we  shall  see  in 
Chapter  VI,  are  well-suited  for  use  in  dynamics. 

DEFINITIONS.  A  dyad  MN  comprises  two  motors  and 
represents  a  linear  transformation  such  that  tor  every  motor 

(MN)  L  =  M(N  o  L) 

n 

A  dyadic  is  a  sum  of  dyads  A  =   J]   MN-  defined  by 

i  =  l 


*   _N  and  J;^*^  would  be  complements,  by  definition,  it  N  n  N^  = 
{0}.  ~   ~ 
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n 
AL  =   y   M.  (N.  o  L) 
i  =  l   '   ' 


We  shall  call  a  dyadic  symmetric  if  it  is  symmetric 
with  respect  to  the  reciprocal  product;  that  is,  for  every 
pair  M,  N,  we  have 

M  o  AN  =  N  o  AM 

(Von  Mises  defined  symmetric  dyadics  likewise.)  We  shall  be 
primarily  concerned  with  symmetric  dyadics,  and  we  now  give 
several  associated  definitions  and  properties  of  symmetric 
dyadics . 

If  the  bilinear  form  M  o  AN  of  a  symmetric  dyadic  A 
is  positive  definite  (i.e.,  M  o  AM  >  0  if  and  only  if 
M  ;«^  0 )  ,  then  A  is  said  to  induce  a  Euclidean  inner  product 

(•|«)   defined  by  (M|N)  =  M  o  an.    The  induced  Euclidean 

1/2 
norm   II  •  n  is   defined   by  IIMII  =  (MJM)  '^  .       All   inner 

products  that  we  shall  consider  are  Euclidean. 

We  call  the  eigenvectors  of  a  dyadic  eigenmotors . 

THEOREM  2.7.   If  a  symmetric  dyadic  A  induces  an  inner 
product  (•!•),   then  there  exists  a  co-reciprocal,  linearly 
independent   set   of   six   eigenmotors  (M  ,  M2 »   ...,   Mg } . 
Moreover,  with  X-,,  X2/   .../  X^  denoting  the  corresponding 
eigenvalues,  the  dyadic  can  be  expressed  in  the  form 


^=  .\    m:^  ''i^    •  (2.3) 

1  =  1   1   1 


26 


PROOF.  Since  A  is  synunetric  with  respect  to  the 
reciprocal  product,  A  is  also  synunetric  with  respect  to 
(•!•);  for,  (M|AN)  =  M  o  (AN)  =  (AM)  o  AN  =  {AM|N).  Thus, 
there  exists  an  orthogonal  set  of  six  eigenmotors  (i.e.,  the 
set  is  orthogonal  with  respect  to  (•!•))/  and  the 
corresponding  eigenvalues  must  be  real.  In  order  to  show 
that   the   set   of   eigenmotors  {M, ,  M. ,  ...,  M^ }  is   co- 

1      Z  D 

reciprocal,  we  note  that 

(M. IM. )  =  M. oAM. 

=  M. oX.M. 
1   D  J 

=  X . (M. oM. ) 
J   1   J 


Now,   since   ('I*)    is   positive   definite,   we   must   have 

Xj  >  0.    Since  ( M^  |  M^  )  =  0  for  i  ^  j  ,  we  have  also  that 

M^  o  M-  =  0. 

In  order  to  prove  the  latter  part  of  the  theorem,  we 

6     X. 
merely  verify  that  the  eigenmotors  of   T  — — rr—  M.M.  are  the 

1  =  1    1   1 

same  as  those  of  A: 


J,     ^H-SW7   "i^)^  =  M-^  M  (M  OM  ) 

1  =  1   11       -•    JJ 


=  X  M. 
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Now,  that  M^  o  M^  7^  0  can  be  proven  two  ways.   Since  {M,  , 
M2 ;   .../   Mg}  is   co-reciprocal   and   linearly   independent, 
by  Theorem   2.6   of   subsection   2.5.2   we   have   that   the 
motors  cannot  be  self-reciprocal.    Alternatively,  M^  o  h^ 

"IT      ^*^i   °    ■^i'^i^  =  X^  ^^i  °  •'^i^  =  IT  ^"il'^i^  >  °-   Q-^-D. 
i  i  i 

6     X. 

We  call  the  expression   T   ..  ^,.   MM.  a  dyad  expansion 

"^  .  ^,  M.  oM.    1  1    — * *^ 

1  =  1   1   1 

for  A.  This  is  similar  to  the  spectral  decomposition  of  a 
linear  transformation;  if  we  were  to  define  the  dyad  M^^M^ 
such    that    (M^Mj^)N   =   Mj^(Mj^|N),    then    the    spectral 

6      X. 

decomposition  of  A  would  be  A  =   V  -rn — m — r  M-M-  . 

i  =  l  ("il'*i)   ^  ^ 

It  is  interesting  to  note  that  if  A  is  symmetric  with 
respect  to  (•!•)/  then  A  always  has  a  spectral 
decomposition;  yet,  if  A  is  symmetric  with  respect  to  the 
reciprocal  product  and  does  not  induce  an  inner  product, 
then  A  might  not  have  a  dyad  expansion.  This  is  because 
eigenmotors  of  A  can  never  be  self-orthogonal;  but  it  is 
possible  that  some  eigenmotors  be  self -reciprocal  (consider 
the  dyadic  MM  where  M  is  self-reciprocal).* 

Theorem  2.7  constitutes  a  generalization  of  some 
results  of  Featherstone  (1984,  see  p.  57).  It  is 
interesting  to  note  that  Featherstone  obtained  the  dyad 
expansion 


*  That  A  induce  an  inner  product  is  only  a  necessary 
condition  for  A  to  have  a  dyad  expansion.  It  is  not 
sufficient  (consider  the  dyadic  MM  where  M  is  an  element  of 
a  co-reciprocal  basis  for  M_,  all  elements  of  which  can  be 
taken  as  eigenmotors  of  MM)  . 
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6    (AM. )(AM. ) 

^         .  S    (M.  M. ) 
1  =  1      1  '  1 


The  dyads,  however,  in  this  dyadic  can  in  fact  be  replaced 
with  the  simpler  expression  in  eqn.  (2.3). 

Finally,  because  of  the  form  of  eqn.  (2.3),  an 
expression  for  the  inverse  A"-'-  is  known  immediately.  Now, 
it  is  easy  to  show  that  A~  induces  an  inner  product  and 
that  the  eigenmotors  of  A  can  also  be  taken  as  the 
eigenmotors  of  A~  .   Thus,  we  have 

_      6     X"-"- 

^"   =  J,  mT^  "i^    •  (2.4) 

1  =  1   1   1 


2 . 6   Motor  Differentiation  and  Integration 

The  existence  of  norms  on  _M  makes  _M_  a  special  kind  of 

metric  space:   a  normed  space.    We  assume  that  limiting 

operations  on  _M  are  with  respect  to  Euclidean  norms,  and  we 

assume   that   limiting   operations  .for   vectors   are   with 

respect  to  the  norm  that  gives  the  magnitude  of  a  vector 

(i.e.,   |v|  =   (v  •  v)-'-/^)^  which  is  also  Euclidean.    We 

consider  motors  as  functions  of  real  variables  only.   Thus, 

with  t,  a  e  _R ,  lim  M(t)   =  L  means  that  for  every  e  >  0, 

t+a 
there  exists  a  5  >  0  such  that  0  <  |t  -  a|  <  6  implies  il  L  - 

M(t)ll    <   e.     Because   normed,   finite-dimensional   vector 
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spaces  are  topoloy ically  closed,*  we  have  that  all  limits 

for  motor  functions  are  motors  (hence,  L  is  a  motor).**   The 

notation  lim  m(  t )  denotes  the  vector  field  such  that  for 
t>a 

every   fixed   point   P   in   the   frame   of   reference,   the 

correspondinu  vector  in  this  vector  field  is  lim  m  (t). 

t>a   i^ 

Before   we   give   the   relationship   between  lim  M(t)    and 

t^a 

lim  m(t),  we  first  need  to  consider  a  few  definitions  and 
t>a 

properties . 

DEFINITIONS.  A  motor  M(t)  is  fixed  in  a  frame  if  for 
every  fixed  point  P  in  this  frame,  mp(t)  is  constant.  A  set 
of  motors  is  fixed  in  a  frame  if  each  motor  in  the  set  is 
fixed . 

We  shall  assume  for  the  remainder  of  this  section  that 
all  motors  are  referred  to  the  same  frame.  Thus,  if  a  motor 
is  fixed,  it  is  fixed  with  respect  to  this  frame. 

Let  {M-  ,  M_  ,  ...,  M^ }  be  a  fixed  basis,  and  let  ai(t), 

a2(t),   ...,  ag(t)  be  the  coordinates  of  M(t).    We  shall 

assume  the  equivalence  of  convergence  in  norm  and  coordinate 

convergence;  that  is,  lim  M{t)  =  0  if  and  only  if  lim  a^lt) 
* t+a  t+  a 


*    A  set  that  contains  its  limit  points  is  said  to  be 
topolog ically  closed  or,  simply,  closed. 


**  See  Lemma  54.2  of  Voyevodin  (1983,  p.  173),  which  states 
that  any  finite-dimensional  subspace  of  a  normed  space  is 
closed . 
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=  0  for  each  i.*   It  follows  that,  on  referring  to  a  fixed 
standard  basis,  we  have 


[lim  M( t) ]  = 
t-»-a 


lim  in(  t ) 

t*a 


lim  mp(t) 
t+a 


We   now   give   the   relationship   between  lim  M(t)  and 

t>a 

lim  m(t).   For  simplicity,  we  suppress  some  notation  so  that 
t  +  a 

lim  M  =  lim  M(t),  lim  m  =  lim  m(t),  etc. 
t+a  t-»-a 

THEOREM   2.7.     For   the   motor   tunction   M(t),   the 

existence  of  lim  M  or  lim  m  implies  the  existence  of  the 

other  and 


lim  M  :  lim  m    , 

for  which  the  principal  vector  is  lim  in. 

PROOF.   Suppose  lim  M  exists,  and  let  M  be  referred  to 
a  standard  basis.   Then, 


[lim  M] 


1  im  m 
lim  m. 


Since  P  is  arbitrary,  we  have 


*  For  a  general  norm,  the  proof  ot  this  is  a  bit 
involved.  For  a  Euclidean  norm  (which  is  all  that  we 
consider),  however,  the  proof  is  straightforward  (Halmos, 
1958,  p.  175). 
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lim  M  :  1 im  m, 

for  which  the  principal  vector  is  lim  fn. 

Suppose  lim  m  exists.   For  every  pair  of  points  P,  Q, 
which  we  assume  are  fixed,  we  have 


mQ  -  rup  =  m  X  PQ 


Taking  the  limit  of  both  sides, 


lim  mQ  -  lim  rap  =  lim(m  x  PQ) 


Since  PQ  is  arbitrary  and  fixed,  we  have 


lim(i5  X  PQ)  =  (lim  m)  x  pQ 


Thus  lim  m  is  a  motor  with  principal  vector  lim  m.   Since 


1  im  m 
lim  mr 


=  [lim  M] 


we  have  lim  m  :  lim  M.   Q.E.D. 

We   use   the   notation  — -  in   order   to   indicate   that 

at 

differentiation   is   taken  with   respect   to   fixed  points; 

.  .     3m    ,  .   m( t ) -m( a ) 
thus  -r—   =    lim  — ^^ — J— — ^ — '-. 

t+a 

COROLLARY.   For  the  motor  function  M(t),  the  existence 

of  either  -r—   H  or  -r?  implies  the  existence  of  the  other  and 

dt       9t    '^ 
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dt  "  •  3t 


for  which  the  principal  vector  is  -3—. 

at 


PROOF.   Apply  the  theorem  to  the  motor  function 
M(t)-M(a)    m(t)-m(a) 


t-a    •     t-a 


.  Q.E.D. 


We  shall  always  assume  that  limits  and  derivatives  of 
motor  functions  exist.  In  addition,  we  shall  assume  that 
motor  functions  are  integrable.  With  these  assumptions,  we 
give  the  following  theorem  without  proof. 

THEOREM  2.8.   For  the  motor  function  M(t),  we  have 

6  S 

/   M(T)dT  :   /   m(T)dT    , 
a  a 

S 
for  which  the  principal  vector  is  /  m(T)dT. 

a 
The  important  results  of  this  section,  as  regards  the 

remainder  of  this  work,  can  be  summarized  as  follows: 


,  ,~      3in_ 

Mt  "'p    Mt'   at' 


6  0       8 

(  /   MdT)p  =  (  /   mdi,  /  fflpdi) 
a  a        a 


2 .7   Differentiation  in  a  Moving  Frame 
The  relationship  between   the   time  derivatives  of  a 

motor  taken  with  respect  to  two  frames  is  analogous  to  that 

d         d ' 
for   a  vector.    If  -^-r ■  M  and  -r-r   M  are  derivatives  taken  with 

at       dt 
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respect  to  Z  and  E'  and  V  the  velocity  motor  of  E'  relative 
to  E,  then 


d       d ' 
dt      dt 


(2.5) 


Before   we   verify   this,   we   need   to   introduce   two 
operations   associated   with   vector   fields:   the   material 


derivative 


dm 


dm 


-nr  r.,    and  the  gradient  7m.    The  notation  -tt-    „, 
dt  E'  -^ dt  E* 

denotes  the  derivative  with  respect  to  E  of  a  vector  m 
associated  with  a  fixed  point  in  E'  (this  vector  field  is 
not  necessarily  a  motor).   The  gradient  7m  at  a  point  P  is  a 
linear  operator  defined  by 
|m--m_,-(  7m)  PQ  i 


lim 
Q+P 


0   P 


=  0 


(2.6) 


PQ 


where  P  can  be  approached  in  any  direction  by  Q.  We  have 
immediately  that  7m  =  mx  (where  we  define  the  operator  inx 
in  the  obvious  way)  since  m  is  a  motor;  for,  with  this 
substitution  in  eqn .  (2.6),  the  numerator  is  zero  for  all 
pairs  of  points  P,0.  Thus  the  gradient  7m  is  independent  of 
the  choice  for  P  and  is  determined  by  the  principal  vector 
m.* 


*   In  vector  analysis,  the  curl  of  a  vector  field  v,  which 
is  denoted  by  7v ,  is  defined  by  the  operator  equation 

(  7v  )  X  =  7v  +  (  7v )  "^ . 

Now,  for  the  motor  M,  (7m)'^  =  -7bi  since  Vm  =  mx.   Thus,  we 
have  that  7xra  =  2n. 
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Now,  by  the  chain  rule  of  differential  calculus,  we 
obtain 


dm 

dt 


dm 


,  =  -g^  +  (  7m)m 


(Note  that  ^  =  ^ 
Vm,  we  have 


.)   On  rearranging  and  substituting  for 


Sml    _  9'm 

3t  E'     3t 


+  V  X  m 


which  is  the  vector  field  form  of  eqn.  (2.5). 


2.8   The  Acceleration  Motor 
The  acceleration  motor  of  a  frame  E'  relative  to  Z    is 
defined  as  the  time  derivative  of  the  velocity  motor  V 
relating  these  two  frames: 


^    dt      3t 


It   is  important  to  note  that  ^  is  not  the  acceleration 

o  t 

of  the   point   with   which   it   is   associated.     Or,   more 
precisely,   it   is   not   the   total   acceleration,   which   is 


dv 


E'* 


The  acceleration 


given  by  the  material  derivative  -r— 
motor  V  comprises  what  is  called  the  local  acceleration 
field  of  E'  relative  to  E.  The  relation  between  the  local 
and  total  acceleration  can  be  determined  from  eqn.  (2.7)  of 
the  previous  section: 
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3t 


dv 

dt 


,  +  V  X  V 


2 . 9   Identities 
Algebraic  Identities: 

M  o  N  =  N  o  M    , 

M  X  N  =  -(M  X  N)     , 

Lo(M+N)=LoM+LoN    , 

Lx(M+N)=LxM+LxN    , 

(XM)  o  N  =  X(M  o  N)     , 

(  XM)  X  N  =  X(M  X  N) 

The  algebraic  identities  can  be  proven  using  algebraic 
properties  of  the  corresponding  vector  fields.   For  example, 

MxN:mxn+nixn 

=    -(nxiii)    -(nxm) 


=    -(n    x    m    +    n    X    in)     :    -(N    x    M) 
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Differential  Identities: 


^    ( XM)     =    XM    +    XM  , 


^    (M   +    N)     =    M   +    N  , 


dt 


(MoN)+MoN    +    MoN  , 


dt 


(MxN)     =MxN+MxN 


The  differential  identities  can  be  proven  using  the 
definition  of  derivative  and  the  algebraic  identities. 
(Note  that  we  need  not  resort  to  the  differential  properties 
of  the  corresponding  vector  fields;  this  is  because  we  have 
defined  limits  tor  motors  with  respect  to  a  metric  on  M_.  ) 
For  example, 


dt 


MxN) 


^  ^^^    M(t)xN(t)-M(a)xN(a) 

t=a    ^  t-a 

t+a 


lim 
t>a 


M(  t)  xN(  t)+M(  t)  xN(a)-M(  t)  xN(  a)-M(  a)  xN(a) 

t-a 


,.    [M( t)-M(a)] xN(a)  ,  ,.   M( t)x[N( t)-N(a) ] 
lim  ; +  iim 


t>a 


t-a 


t*a 


t-a 


=  MxN+MxN 
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The  following  additional  identity  proves  very  useful  in 

rigid  body  dynamics  (see  Chapter  VI): 

LoMxN=LxMoN 


by 


CHAPTER  III 
GEOMETRICAL  CONSIDERATIONS 


3.1   The  Pitch  and  Central  Axis 


DEFINITION.   The  pitch  of  a  non-zero  motor  M  is  given 


r 


ni'm^ 


ni*m 


h  = 


if 


if 


m  ^  0 


m  =  0 


A  non-zero  motor  is  said  to  be  proper  if  h  is  finite 
and  improper  if  h  is  infinite.  The  product  in  *  mp  is 
independent  of  the  choice  of  point  P  as  can  be  verified  from 
n.  (2.1).   In  fact,  fn  *  mp  =  ^  M  o  m   (for  all  motors);  so 


eqr 


we  can  alternatively  write  h  =  ^  ~z~Z   ^^'^  proper  motors. 

m  •m 
THEOREM  3.1.    For  every  proper  motor,  there  exists  a 

unique  line  such  that  rap  =  hm  for  every  point  P  on  the  line. 

PROOF.    We  assert  that  the  equation  of  the  line  such 

that  rap  =  hm,  P  being  a  point  on  the  line,  is  given  by 


mxn 


r  = 


o 


+  km 


k  e  R 


m< 
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where  0  is  an  arbitrary  point  assumed  not  on  the  line  ana 
r  =  OP. 

By  definition  of  a  motor,  we  have 

fflp  =  mQ  +  m  X  r 

Substituting  for  r  and  simplifying. 


mxm     ^ 

mr,  =  m^  +  m  X  ( +  km) 


=  hIq  +  ( hm  -  mQ  +  0 ) 


=  hm 


mxm 

which  verifies  the  assertion.    Now,  r  =  for  k  =  0. 

m  'm 
Since   points   O  and   P  are  distinct   (by  assumption),   we 

mxm 

have    ^  ^   *    0,   which  in  turn  implies  m    x   m^   i^     0.    Thus, 

m*m 
mQ  7^  hm  (for,  mQ  =  hm  implies  in  x  mQ  =  0).   Thus  the  given 

line  is  unique.   Q.E.D. 

The  line  for  a  proper  motor  such  that  rap  =  hm  for 
points  P  on  the  line  is  called  the  central  axis  of  the 
motor . 

THEOREM  3.2.  Let  M  be  proper.  Then,  every  non-zero 
motor  in  <M>  has  the  same  central  axis  and  pitch.  Moreover, 
only  motors  in  <M>  have  this  central  axis  and  pitch 
combination . 
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PROOF.   Let  P  be  on  the  central  axis  of  M.   Then, 


mp  =  hm  (3.1) 


where  h  is  the  pitch  of  M.   Multiplying  by  a  non-zero  scalar 


Xntp  =  hXm    / 


and  noting  that 

(XM)p  =  {Xm,     Xmp)    , 

we  conclude  that  P  is  also  on  the  axis  of  XM  and  that  the 
pitch  of  XM  is  also  h. 

For  the  latter  part  of  the  theorem,  suppose  motor  N  has 
the  same  central  axis  and  pitch  as  M.   Then, 


Hp  =  hn    .  (3.2) 


Since  n  and  fn  are  parallel,  there  exists  a  scalar  3  such 
that  n  =  em.  Thus,  from  eqns.  (3.1)  and  (3.2),  we  also  have 
Hp  =  Sfflp.   Now, 


( N) p  =  (n,  np) 
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=  (Sin,  gmp) 

=  {eM)p 

Therefore,  N  =  BM  e  <M>.   Q.E.D. 

Thus,  we  have  that  one-dimensional  subspaces  of  proper 
motors  characterize  central  axis  and  pitch  combinations. 

For  improper  motors,  there  is  no  central  axis  in  the 
sense  defined  above.  We  note  that  an  improper  motor  M 
comprises  a  constant  vector  field  (i.e.,  for  every  pair  of 
points  P,  0,  nip  =  hIq)  .  Since  in  =  O  and  h  =  <»,  there  can  be 
no  points  P  such  that  mp  =  hm.  We  now  resort  to  projective 
geometry  in  order  to  define  the  central  axis  of  an  improper 
motor. 

Since  one-dimensional  subspaces  of  proper  motors 
characterize  central  axis  and  pitch  combinations,  we 
likewise  define  the  central  axis  of  an  improper  motors 
(which  by  definition  has  infinite  pitch)  as  some  geometric 
entity  characterized  by  the  one-dimensional  subspace  to 
which  the  motor  belongs.  In  projective  geometry,  this  new 
type  of  central  axis  is  known  as  a  line  at  infinity  or  a 
non-Euclidean  line,  as  opposed  to  a  Euclidean  line,  which  is 
a  line  in  Euclidean  space. 

By  Theorem  3.2  and  the  above  definition  of  the  central 
axis  for  a  improper  motor,  we  can  state  the  following: 
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THEOREM  3.3.  There  exists  a  one-one  correspondence 
between  one-dimensional  subspaces  of  M_  and  all  central  axis 
and  pitch  combinations. 

3 .2   On  Representing  Lines,  Rotors  and  "Vectors" 

By  definition,  we  have  that  for  every  non-Euclidean 
line,  there  exists  a  motor  (which  is  improper)  that  has  this 
line  as  a  central  axis.  We  must,  however,  establish  an 
analogous  fact  regarding  Euclidean  lines. 

For  some  arbitrary  Euclidean  line,  let  v  be  a  non-zero 
vector  bound  to  the  line  and  r  be  directed  from  an  arbitrary 
point  in  space  to  any  point  on  the  line.  Then,  the  moment 
distribution  r  x  v  is  a  motor  with  principal  vector  v.  (If 
V  is  a  force,  then  this  motor  is  a  wrench.)  Let  M:r  x  v. 
Then,  in  =  v  and  m  =  r  x  v.  The  pitch  h  of  M  is  clearly  zero 
since  m  •  m  =  0.  Thus,  we  have  that  the  given  line  is  the 
central  axis  of  M  since  for  any  point  P  on  the  line  mp  =  hm 
(=0),  and  we  can  now  state  the  following: 

THEOREM  3.4.  For  every  line  (Euclidean  and  non- 
Euclidean),  there  exists  some  motor  that  has  this  line  as  a 
central  axis. 

We  can  thus  represent  every  line  (Euclidean  and  non- 
Euclidean)  by  a  motor.  Therefore,  we  can,  in  turn, 
represent  a  line  by  motor  coordinates  or  by  a  bivector. 
When  a  motor  M  is  of  zero  or  infinite  pitch,  then  the 
coordinates  of  M  with  respect  to  a  standard  basis,  that  is. 
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[M]  = 


m 


constitute  the  well-known  Plucker  coordinates  of  the  line* 
that  is  the  central  axis  of  M. 

DEFINTITION.   A  rotor  is  a  motor  of  zero  pitch. 

The  term  rotor  is  adopted  from  Clifford  (1873),  which 
he  defined  in  the  sense  of  a  line  vector  (i.e.,  a  vector 
bound  to  a  Euclidean  line),  and  the  definition  we  use  is 
equivalent  to  Clifford's.  The  discussion  preceding 
Theorem  3.4  shows  that  there  is  a  unique  motor  of  zero  pitch 
with  a  given  principal  vector  and  central  axis  (here,  we  can 
consider  the  principal  vector  bound  to  the  central  axis), 
and  it  turns  out  that  our  definition  for  addition  of  rotors 
(which,  of  course,  is  simply  addition  of  motors  of  zero 
pitch)  is  also  equivalent  to  Clifford's. 

Rotors  can  thus  be  used  to  represent  any  Euclidean  line 
and/or  line  vector.  (In  particular,  rotors  can  be  used  to 
represent  forces  when  forces  are  interpreted  as  line 
vectors  .  ) 

Now,  what  we  call  an  improper  motor  is  what  Clifford 
defined  a  vector  to  be.   Thus,  a  "vector"  is  the  set  of  all 


*   The  Plucker  coordinates  of  a  line  are  called  homogeneous 
coordinates.    Coordinates  of  a  geometric  entity,  say,  a-,. 


,  a^      are   homogeneous   when   Xai,  Xa-y 


represent  the  same  geometric  entity  where"  X  is  "an  arbitrary 
non-zero  real  scalar. 
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vectors  indexed  to  points  in  space  that  have  the  same 
magnitude  and  direction  (i.e.,  a  constant  vector  field). 
Addition  of  "vectors,"  as  Clifford  defined  it,  is  equivalent 
to  addition  of  improper  motors. 

As  Clifford  pointed  out,  the  set  of  "vectors"  is  closed 
under  addition  whereas  the  set  of  rotors  is  not.*  Clifford 
designated  a  motor  as  that  which  may  be  a  rotor,  "vector," 
or  the  new  entity  that  can  result  from  the  addition  of  two 
rotors.  This  new  entity  is  what  we  call  a  motor  of  non- 
zero, non-infinite  pitch.  We  shall  show  in  the  next  section 
that  such  a  motor  can  be  decomposed  into  a  rotor  and  an 
improper  motor  (i.e.,  a  "vector"). 

We  close  this  section  with  a  geometric  description  of  a 
standard  basis.  Let  {E,  ,  E_  ,  ...,  E, }  be  such  a  basis  for 
which 


(E^)q  =  (e^,  0),  (£2)0  =  (^2'    °^'  ^^3^0  "  ^®3'  °^ 


(£4)^  =  (0,  e-,^),  (£5)^  =  (0,  62),  (Eg)y  =  (0,  03) 

We  note  that  E-,,   E2  /   E-,  are  rotors  and  Ea,      Ec  ,  Ec  are 
improper  motors.   The  central  axes  of  Ei,  E2,  E-j  constitute 


*  We  tacitly  assume  here  that  the  zero  "vector"  and  zero 
rotor  are  defined;  these  would  both  be  the  same  as  the  zero 
motor.  Note,  however,  that  we  shall  maintain  the  convention 
that  the  zero  motor  is  neither  proper  nor  improper  (or, 
equivalently ,  the  pitch  is  undefined). 
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a  rectangular  set  of  lines  that  intersect  at  0  (see 
Fig.  3.1).  If  we  designate  the  directions  of  these  lines  by 
e-,  ,  e2i  e^/  we  then  have  a  Cartesian  coordinate  system  with 
origin  0.  Clearly,  the  first  three  coordinates  of  a  motor  M 
with  respect  to  this  standard  basis  are  the  vector 
coordinates  of  m  in  this  Cartesian  system  and  the  last  three 
coordinates  are  the  vector  coordinates  of  m,-.  in  this  same 
Cartesian  system. 

3 . 3   A  Useful  Decomposition 
Let  M  be  proper  with  pitch  h.   We  define 


M^  :m  -  hm     and      M°°:m 


so  that  we  have 


M  =  M°  +  hM"    .  (3.3) 


The   superscripts   indicate   the   pitches.  We   can   show 

simultaneously  that  the  pitch  of  M^  is  in  fact  zero  and  that 

the  central  axis  of  M^  is  the  same  as  that  of  M:   It  P  is  a 
point  on  the  central  axis  of  M,  then 


0  «,     ~ 

nip  =  fflp  -  hm  =  Om 


where  m^  :  M*^  . 


46 


Fig.  3.1, 
Basis 


The  Cartesian  Coordinate  System  of  a  Standard 
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We  thus  have  that  every  proper  motor  can  be  decomposed 
into  a  rotor  and  an  improper  motor  {i.e.,  a  "vector"  in  the 
sense  of  Clifford  as  mentioned  in  the  preceeding  section). 

Now,  from  the  previous  section  we  have  that 


M°:r  X  m 


where  r  is  directed  from  an  arbitrary  point  in  space  to  the 
central  axis  of  m"^  .   Thus,  we  have 


M :  r  X  in  +  hin    ,  (3.4) 

which  enables  us  to  construct  a  proper  motor  given  the 
central  axis,  pitch,  and  principal  vector. 

Finally,  we  note  that  with  respect  to  a  standard  basis, 
[M^]  comprises  the  Plucker  coordinates  of  the  central  axis 
of  M. 

3 . 4   A  Note  on  the  Dual  Operator  and  Dual  Numbers 
It  is  perhaps  of  interest  here  to  consider  a  special 

operator  for  motors  introduced  by  Clifford  (1873),  which  he 

designated  by  o)  so  defined  that 

u)M  =  m"    and    ojCojM)  =  0 
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for  every  proper  motor  M.  It  is  easy  to  show  that  the  rank 
of  the  linear  transformation  represented  by  to  is  three;  the 
null  space  comprises  all  improper  motors  (and  the  zero 
motor).  (That  ojN  =  0  for  improper  N  follows  from  the  fact 
that  N  =  a)((n,  v))"-*-  where  v  can  be  any  vector  and  P  any 
point.  The  motor  ((n,  v))p-'-  is,  of  course,  necessarily 
proper  since  n  ^  0.)  With  respect  to  a  standard  basis,  the 
matrix  of  this  transformation  is 


0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

1 

0 

0 

0 

0 

0 

0 

1 

0 

0 

0 

0 

0 

0 

1 

0 

0 

0 

We  can  equivalently  define   oi  as  a  special  kind  of 
number  such  that  for  all  a,  B  g  R  and  every  proper  motor  M, 


(a  +  6u  +  w^)M  =  ctM  +  6M 


(3.5) 


From  this,  we  can  conclude  oj^  =  0  (chose  a  =  3  =  0)  and  uO  = 
0  (choose  6  =  0).  (Note  that  since  cj^M  =  (o(  wM)  ,  we  have 
that  (jjN  =  0  for  any  improper  motor  N.  )  Thus,  multiplication 
of  a  motor  by  the  number  o)  is  equivalent  to  operating  on  it 
by  the  operator  u. 
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Numbers  o£  the  torm  a  +  Bai  were  caliea  dual  numoers  by 
Stuay  (1903)  and  were  applied  by  him  in  linear  geometry. 
The  term  dual  operator  is  otten  used  in  recent  literature 
for  cu  when  the  operator  interpretation  is  used. 

The  alyebra  of  dual  numbers  as  induced  from  e^n.  (3.5) 
is  a  commutative  riny .  Some  authors  who  have  investigated 
properties  of  dual  number  alyebra  include  Dimentbery  (i96b), 
Veldkamp  (1976),  and  Rooney  (1975).  The  motivation  tor  the 
study  of  dual  number  alyebra  lies  primarily  in  its 
application  to  line  geometry  and  rigid  body  kinematics 
(Study,  1903;  Brand,  1947;  Dimentbery,  1965;  Yany ,  1969, 
1974;  Yany  and  t'reudenstein ,  1964;  Veldkamp,  1976;  Keler, 
1973,  1979).  The  application  of  dual  number  algebra  to 
rigid  body  dynamics,  however,  has  received  much  less 
attention  (Dimentbery,  196b;  Yany,  1967,  1969,  1971,  1974; 
Pennock  and  Yany,  1983). 

Dual  number  algebra  is  not  a  particularly  useful  tool 
for  rigid  body  dynamics.  Although  equations  ot  motion  can 
be  expressed  in  dual  number  torm  (see  tne  cited  works  ot 
Dimentberg  and  Yang),  such  a  form  does  not,  in  my  opinion, 
provide  a  catalyst  tor  explication  ot  dynamic  properties. 
It  is  for  this  reason  that  we  will  not  exploit  dual  number 
algebra . 
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3  . 5  Streamlines  of  a  Motor 
By  use  of  the  decomposition  in  eqn.  (3.3),  we  can 
easily  describe  the  geometric  nature  of  the  vector  field  of 
a  motor.  Let  M  be  proper  and  let  r  =  OP  where  0  is  any 
point  and  P  is  on  the  central  axis  of  M.  Then  we  have  mx  = 
r  X  m  (where  m^:M'^).  The  streamlines  of  m^  comprise 
concentric  circles  about  the  central  axis  and  in  planes 
normal  to  the  central  axis.  One  such  circle  is  illustrated 
in  Fig.  3.2.  By  adding  hm  to  m^  (note  hm  +  m  :M),  the 
streamlines  for  m  become  concentric  helices  of  pitch  h  about 
the  central  axis.  One  such  helix  is  illustrated  in 
Fig.  3.3.  The  streamlines  of  m°°(  =m)  are  simply  straight 
lines  parallel  to  the  central  axis. 

3.6   The  Unit  Motor 

We  reserve  the  symbol  S  for  unit  motors.  A  unit  motor 
S  is  defined  such  that  \s\  =  1  (|s|  =  1)  if  S  is  proper 
( improper ) . 

It  is  important  to  note  that  a  unit  motor  is  not 
defined  with  respect  to  a  norm.  The  reason  for  so  defining 
unit  motors  is  that  they  are  convenient  for  geometric 
purposes  and  for  use  in  mechanics.  (Von  Mises  defined  unit 
motors  similarly.) 
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r    X   m 
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Fig.    3.2.    Streamline   of   a   Rotor 


^.    r    X   in  +    lim 


m 


Fiq.    3.3.    Streamline    of    a   Motor    (hT^O    and    h?^") 
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3.7   The  Screw 

DEFINITION  (Synthetic).  A  screw  is  a  central  axis  and 
pitch  combination  of  some  motor. 

DEFINITION  (Analytic).  A  screw  is  a  one-dimensional 
subspace  of  M_« 

The  terms  synthetic  and  analytic  are  taken  from 
projective  geometry.  The  synthetic  definition  of  the  screw 
is  the  same  as  that  given  by  Ball  (1900).  We  shall, 
however,  usually  confound  the  two  definitions  by  exploiting 
Theorem  3.3.  The  synthetic  definition  is  convenient  for 
visualization  purposes,  and  the  analytic  definition  is 
convenient  for  algebraic  purposes. 

In  analogy  with  motors,  we  say  that  a  screw  is  proper 
( improper)  if  its  pitch  is  finite  (infinite). 

We  can  now  define  a  proper  motor  in  terms  of  a  screw. 
It  is  in  this  manner  that  Clifford  defined  the  motor:  "Just 
as  a  vector  ...  is  magnitude  associated  with  direction, 
and  as  a  rotor  ...  is  magnitude  associated  with  an  axis; 
so  this  new  quantity,  which  is  the  sum  of  two  or  more  rotors 
.  .  .  ,  is  magnitude  associated  with  a  screw."  It  is  in 
fact  implicit  in  Clifford's  definition  of  motor  (and  rotor) 
that  a  sense  is  assumed  for  the  central  axis.  Thus,  the 
central  axis,  pitch,  and  the  magnitude  and  sense  (with 
respect  to  the  central  axis)  of  the  principal  vector 
determines  a  proper  motor. 
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We  can  now  give  Ball's  definitions  of  twist  and 
wrench.  A  "twist  about  a  screw"  is  equivalent  to  a  velocity 
motor.  The  twist  is  the  magnitude  of  the  angular  velocity 
(or  linear  velocity  if  the  velocity  motor  is  improper).  A 
"wrench  on  a  screw"  is  the  same  as  our  definition  of 
wrench.  Ball's  wrench  is  the  magnitude  of  the  net  force  of 
the  force  system  (or  the  couple  of  the  force  system  if  the 
representing  wrench  is  improper).  (As  with  Clifford,  a 
sense  is  assumed  for  the  central  axes  of  screws.) 

By  the  use  of  unit  motors,  we  can  make  the  above  notion 
of  "magnitude  associated  with  a  screw"  precise.  Every  motor 
can  be  expressed  as  a  scalar  multiple  of  a  unit  motor: 
M  =  XS  where  X  =  ±|fn|  (X  =  ±|in|  for  improper  M)  and  S  =  5! 

A 

(S  can  be  chosen  as  any  unit  motor  if  X  =  0).  Thus,  S 
represents  the  screw  of  M  together  with  a  direction,  where 
we  assume  the  direction  of  S  is  given  by  s  (s)  if  S  is 
proper  (improper).  We  assume  there  is  no  definite  direction 
if  X  =  0.  We  call  X  the  magnitude  of  the  motor,  which, 
obviously,  may  be  positive  or  negative  (we  note  in 
particular  that  the  magnitude  of  each  unit  motor  S  may  be 
either  ±1).  If  M  is  a  wrench,  then  X  is  the  magnitude  of  a 
force  (torque)  if  M  is  proper  (improper);  if  M  is  a  velocity 
motor,  then  X  is  the  magnitude  of  an  angular  (linear) 
velocity  if  M  is  proper  (improper). 

The  following  theorem  is  useful  in  kinematics. 
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THEOREM  3.5.  A  unit  motor  S  is  fixed  in  a  frame  E  if 
and  only  if 

(i)      the  central  axis  of  S  is  fixed  in  E; 
(ii)     the  pitch  of  S  is  constant; 
(iii)    the  direction  of  S  is  constant. 

PROOF.  The  proof  is  immediate  for  improper  S.  If  S  is 
proper,  then  conditions  (i)  and  (iii)  determine  S^,  and  S*^ , 
which  in  turn,  determines  s".  Both  S  and  S"  must  be 
fixed.  That  the  pitch  h  is  constant  implies  that  S  =  S  + 
hs"  is  fixed. 

The  proof  of  the  converse  is  immediate.   Q.E.D. 

DEFINITIONS.  The  screws  <M>  ,  <N>  are  said  to  be 
reciprocal  to  each  other  (or,  simply,  reciprocal)  if  M  o  n  = 
0  (i.e.,  the  motors  M,  N  are  reciprocal).  A  motor  system  is 
a  subspace  of  _M.  A  screw  system  is  the  set  of  screws 
corresponding  to  a  motor  system.  The  order  of  a  motor 
system  and  its  corresponding  screw  system  is  the  dimension 
of  the  motor  system.  A  one-,  two-,  ...,  n-system  is  a  screw 
system  of  order  one,  two,  ...,  n. 

The  term  screw  system  is  due  to  Ball  and  the  term  motor 
system  is  due  to  Everett  (1875). 

Motor  systems  are  convenient  for  characterizing 
permissible  relative  motion  between  rigid  bodies  in 
constrained  rigid  body  systems  (see  Chapter  V). 
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3 .8  Geometry  of  Motor  Operations 
In  this  section,  we  present  some  geometric  properties 
associated  with  the  reciprocal  product,  cross  product,  sum, 
and  linear  combination.  All  the  results  obtained  for  the 
first  three  operations  are  well-known.  The  derivation  of 
the  circle  representation  for  a  two-system  in  subsection 
3.8.4,  however,  is  believed  to  be  simpler  than  any  other 
derivations  found  in  the  literature. 

In  the  following,  we  assume  that  the  pitches  of  the 
motors  M-,,  M2 ,  and  M  (and  the  unit  motors  S-,  ,  82^  and  S)  are 
h-i  ,  hp ,  and  h . 

3.8.1   The  Reciprocal  Product 

Let  Ml  =  ^iS-1  and  M2  =  ^9^2  ^^  proper,  the  central  axes 
of  which  are  illustrated  in  Fig.  3.4.  The  common  normal 
line  intersects  the  central  axes  at  points  P-,  and  P2.  We 
specify  the  relative  position  of  one  axis  relative  to  the 


other  by  the  vector  de  =  P,P^,      where  |e|  =1  and  d  >    0,    and 
by  the  angle  a,     which  is  subtended  by  s-,  and  S2  and  is 
measured  in  a  right  hand  sense  about  e. 
Consider 


M-|  o  M2  =  X,X2St  o  S2 


=  X^X2(S°   +  h^s")  o  (3°   +  h2S2) 
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Fig.    3.4.    Central   Axes   of   S,    and   S- 
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X3^X2(S^    ^    S^    +    h^S^    °    ^2    +    h^S^    o    sj^    +    ^1^2^2    °    ^2       ' 


Now, 


sJosO 


'=?'Pj    °    '^"'P^ 


(s^,    0)    o    (S2,    de    X   s^) 


=   -de    •    s,    ^    s^ 


=    -de    •    esina 
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=   -d    sin a 
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00  0         «.  ^ 

^1    °    ^2    =    ^1    •    ^2 


=    cosa 
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OO  00  -V  .N. 


=    cosa 


12 


ea  OB 

S^    o    S^    =    0  . 


58 


Substituting,    we   obtain 


M^    o   M      =    X    X    [(h      +   h    )cosot   -    d    sina]  ,  (3.6) 


which   is   the   well-known   expression   for   the  reciprocal 

product  of  two  proper  motors.* 

Using  eqn.  (3.6),  we  can  deduce  the  following: 

THEOREM  3.6.   Suppose  the  proper  motors  M-i  ,  M2  (screws 

<M-|^>,  <M2>)    have  intersecting  central  axes.  Then,  the 

motors  (screws)  are  reciprocal  if  and  only  if  one  or  both  of 

the  following  hold: 


(i)   h-j^  =  -h2, 

(ii)  the  central  axes  are  perpendicular. 


It  is  clear  that  every  pair  of  improper  motors  (screws) 
is  reciprocal.  If  M-j^  is  proper  and  M2  improper,  then  it  is 
also  clear  that  M-j^  o  M2  =  0  if  and  only  if  m-^    •  m2  =  0. 

3.8.2   The  Cross  Product 


THEOREM  3.7.    Let  M^ ,  M2  ^^  proper  motors  with  non- 
parallel  central  axes.    Then,  M,  '^  ^2  is  a  proper  motor 


*   Ball  called  the  expression 

— [(hj^  +  h2)  cosa  -  d  sina] 

the  virtual  coefficient  of  the  screws  <S,>  and  <S2>. 
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whose  central  axis  is  the  common  normal  of  the  central  axes 
of  Ml  and  M2 . 

PROOF.    The  product  M-,  x  M2  is  proper  since  m^    x  ^2 
^    0.   Now,  since 


and 


in-]  X   m2  *  Si "  =  m-i  x  ^2  *  ni2  =  0 


M-,  X  M2  o  M-j^  =  M-j^  X  M2  o  M2  =  0 


we  have  by  Theorem  3.6  that  the  central  axis  of  M-j^  x  M2  is 
the  common  normal  to  the  central  axes  of  M-j^  and  M2 .   Q.E.D. 

THEOREM  3.8.  The  product  M-j^  x  M2  is  improper  if  and 
only  if  either  condition  holds: 

(i)   Ml  is  improper  and  m-,  x  £2  /^  0, 

(ii)  Ml  and  M2  are  proper  motors  whose  central  axes 
are  distinct  and  parallel. 

PROOF.   The  proof  is  immediate.   Q.E.D. 

3.8.3   The  Motor  Sum 

THEOREM  3.9.  If  the  proper  motors  M^^  ,  M2  have  non- 
parallel  central  axes  (or  parallel  central  axes  such  that 
iSi  ^  -£2)  ,  then  the  sum  M  =  M-j^  +  M2  is  a  proper  motor  whose 
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central  axis  intersects  the  common  normal  line  (or  every 
common  normal  line)  of  the  central  axes  of  M-.  and  M2 . 

PROOF.  Let  the  common  normal  line  (or  some  common 
normal  line)  be  represented  by  the  rotor  S.  Then,  we  must 
have 

s  *  £■]  =  s  •  m2  =  0    , 

and   by   Theorem    3.6,    we    have 


SoM-|^    =    SoM2=0 


Now,  since  m-,  +  ra2  7^  0  (thus  M  is  proper)  ,  we  have  that  the 
central  axis  of  M  is  perpendicular  to  central  axis  of  S; 
for, 

s    (m-j^  +  m2)  =  s  *  m-j^  +  s    m2 


=  0 


In  addition,  since 


S  o  (M^  +  M2)  =  S  o  M]^  +  S  o  M^ 


=  0 
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it  follows  from  Theorem  3.6  that  the  central  axis  of  M  also 
intersects  the  central  axis  of  S.   Q.E.D. 

THEOREM  3.10.  Let  M-|^ ,  M2  be  proper  motors  whose 
central  axes  are  parallel  and  m-,  =  -ni2 .  Then,  the  sum  M  = 
M-j^  +  M2  is 

(i)  improper  if  the  central  axes  of  M-i  and  M2  are 
distinct  or  h-j^  ^  h2/  or,  equivalently ,  <M-i>  ^ 
<M2>; 

(ii)  zero  if  the  central  axes  of  M-,  and  M2  are  the 
same  and  h-j^  =  h2/  or,  equivalently,  <Mn>  =  <M2>. 

PROOF.  For  both  (i)  and  (ii)  we  have  m  =  m-,  +  m2  =  0. 
For  (i),  <M-|^>  ^  <M2>  implies  that  M  =  M-j^  +  M2  7^  0;  thus,  M 
is  improper.  For  (ii),  <M-|^>  =  <M2>  and  m-,  =  -^a.^  implies 
that  M^  +  M2  =  0. 

THEOREM  3.11.  If  two  proper  motors  W^,  Vi^  have  non- 
parallel,  coplanar  central  axes,  then  the  central  axis  of 
the  sum  M  =  M-j^  +  M2  passes  through  the  point  of  intersection 
of  the  central  axes  of  M-j^  and  M2  if  and  only  if  h-,  =  h2. 
Moreover,  h  =  \\^{.=  h2 )  where  h  is  the  pitch  of  M,  and  the 
central  axis  of  M  lies  in  the  plane  of  the  central  axes  of 
M-,  and  M2 . 

PROOF.  Let  P  be  the  point  of  intersection  of  the 
central  axes  of  M^^  and  M2 .   Then,  we  have 
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(M-|^  +  M2)p  =  (M;i^)p  +  (M2)p 


=  (m-j^,  h-j^m-]^)  +  (£2^  i^2™2^ 

=  (m-i  +  va.2 1    h-iin-i  +  h-jiii2)    • 

Now,  M  +  N  is  proper  since  m-,  +  m2  7^  0.  Point  P  is  on  the 
central  axis  of  M  if  and  only  if  (m-^  +  £2)  ><  (h-,  m-^  +  '^2™2^ 
=  0,  in  which  case  h-,  =  h2  so  that  Si  +  £2  =  h-,(m-|  +  m-i  ) 
(=  h2(i5-i  +  ra2^  ^  ^"^  ^  ~  h-,(=  h2). 

That  the  three  central  axes  are  coplanar  follows  from 
Theorem  3.9.   Q.E.D. 

3.8.4   The  Linear  Combination 


In  this  subsection,  we  study  the  geometry  of  the  linear 
combination  of  two  motors;  that  is,  we  study  the  geometry  of 
the  two-system.  (We  are  assuming,  of  course,  that  the  two 
motors  constitute  a  linearly  independent  set.) 

THEOREM  3.12.  Every  pair  of  distinct  screws  in  a  two- 
system  determines  the  two-system. 

PROOF.  Let  <S-,>  and  <S2>  be  distinct  screws  in  a  two- 
system.  Then,  <S-|>  n  <S2>  =  {0}  implies  that  {S  ,  S  }  is 
linearly  independent  and  thus  forms  a  basis  for  the  motor 
system  that  corresponds  to  the  two-system.   Q.E.D. 

(Note  that  this  theorem  cannot  be  extended  to  screw- 
systems  of  order  greater  than  two:  it  is  possible  that  <Sj^>, 
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<S2>/  •../  <Sj^>  be  distinct  and  that  {S,  ,  S  ,  ...,  S  }  is 
linearly  dependent  for  n  >  2.) 

DEFINITION.  The  ruled  surface  comprising  the  central 
axes  of  a  two-system  is  called  a  cylindroid . 

The  cylindroid  was  discovered  by  Hamilton  (1830)  and 
named  by  Cayley.* 

A  two-system  is  characterized  by  its  cylindroid  and  the 
distribution  of  pitch  for  the  lines  on  the  cylindroid. 

We  classify  two-systems  according  to  the  following 
scheme : 

CLASS  1.   All  screws  are  proper. 
CLASS  2.   All  screws  are  improper. 
CLASS  3.   Only  one  screw  is  improper. 

THEOREM  3.13.  The  above  classification  scheme  is 
exhaustive  and  each  class  is  non-empty. 

PROOF.  Let  <S-,>  and  <S2>  be  distinct  screws  in  a  two- 
system.  That  this  classification  scheme  is  exhaustive 
follows  from  the  fact  that  if  <S-.>  and  <S2>  are  improper, 
then  all  screws  in  the  two-system  are  improper. 


*  In  the  footnote  on  p.  20  of  Ball's  work,  he  writes,  "The 
name  cylindroid  was  suggested  by  Professor  Cayley  in  1871  in 
reply  to  a  request  I  made  when  in  ignorance  of  the  previous 

work  of  both  Plucker  and  Battaglini,  I  began  to  study  this 
surface . " 
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We  must  show  that  each  class  is  non-empty.  Clearly, 
there  exist  two-systems  that  belong  to  Class  2.  Now, 
suppose  <S-|>  and  <S2>  are  proper.  If  the  central  axes  of 
these  screws  are  non-parallel,  then  <S-,>  and  <S2>  determine 
a  two-system  that  comprises  proper  screws  (i.e.,  a  Class  1 
two-system)  .  If  the  central  axes  of  <S-^>  and  <S2>  are 
parallel  and  s-,  =  -S2 ,  then  the  one  improper  screw  of  the 
Class  2  two-system  determined  <S-^>  and  <S2>  is  <S-^  +  S2>. 
Q.E.D. 

We  note  that  every  pair  of  distinct  screws  in  a  Class  1 
two-system  must  have  non-parallel  central  axes;  for,  a  two- 
system  that  has  a  pair  of  distinct  screws  whose  central  axes 
are  parallel  must  also  have  an  improper  screw  (see  Theorem 
3.10).  We  also  note  that  the  proper  screws  of  a  Class  3 
two-system  must  have  parallel  central  axes:  the  central  axes 
of  a  proper  motor  and  the  sum  of  this  motor  and  any  improper 
motor  are  parallel.  By  Theorem  3.9,  we  also  have  that  the 
central  axes  of  the  proper  screws  of  a  Class  3  two-system 
are  coplanar. 

DEFINITION.  A  two-system  is  called  degenerate  if  it 
contains  at  least  one  improper  screw. 

Thus,  Class  2  and  Class  3  two-systems  are  degenerate, 
and  Class  1  two-systems  are  non-degenerate. 

The  cylindroids  of  degenerate  two-systems  are  simple 
geometrically.   The  central  axes  of  a  Class  2  two-system  all 
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lie  in  the  so-called  plane  at  infinity.*  The  central  axes 
of  the  proper  screws  of  a  Class  3  two-system  are  either  all 
the  same  or  lie  in  a  plane. 

The  cylindroids  of  non-degenerate  two-systems,  however, 
are  much  more  complex,  and  we  devote  the  remainder  of  the 
subsection  to  non-degenerate  two-systems.  We  shall 
ultimately  obtain  the  so-called  circle  representation  of  a 
non-degenerate  two-system,  which  was  discovered  by  Lewis 
(1880)  and  used  to  great  advantage  by  Ball  (1900).  The 
circle  representation  quite  elegantly  illustrates  the 
geometric  properties  of  non-degenerate  two-systems  as  can  be 
evidenced  in  Ball's  work.  Since  the  circle  representation 
applies  only  to  non-degenerate  two-systems  (this  is  the 
motivation  for  the  above  definition),  we  shall  assume 
henceforth  that  "two-system"  means  "non-degenerate  two- 
system"  and  that  cylindroids  are  of  two-systems  (i.e.,  non- 
degenerate  ones). 

THEOREM  3.14.  All  lines  on  a  cylindroid  intersect  a 
common  normal  line. 

PROOF.  Every  pair  of  lines  on  a  cylindroid  has  a 
unique  common  normal  line.  It  follows  from  Theorem  3.9  that 
there  exists  a  unique  common  normal  line  for  all  lines  on 
the  cylindroid.   Q.E.D. 


*    In  projective  geometry  all  lines  at  infinity  and  all 
points  at  infinity  lie  in  the  plane  at  infinity. 
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The  common  normal  line  of  a  cylindroid  is  called  the 
central  axis  of  the  cylindroid. 

LEMMA.  To  each  line  on  a  cylindroid  there  exists  a 
corresponding  unique  line  on  the  cylindroid  that  is 
perpendicular  to  it. 

PROOF.  Let  <S-]_>  and  <S2>  be  distinct  screws  of  a 
two-system.  A  screw  in  this  two-system  whose  central 
axis  is  perpendicular  to  the  central  axis  of  <S-i>  is 
<S2  -(s  •   S2)  S-^> .  This  perpendicular  central   axis   is 

unique;  for,  otherwise,  there  would  be  two  parallel  central 
axes,  which  is  impossible.   Q.E.D. 

If  <S-,>  and  <S2>  are  distinct  screws  of  a  two-system, 
every  screw  in  this  two-system  can  be  given  by  <c  .S-,  +  s  ,S-)> 
where  c,  =  cosi|)  and  s.  =  simj)  and  for  which  ipe[0,  tt).  (This 
follows  since  it  is  the  ratio  s^:c.  that  determines  the 
screw,  and  all  possible  ratios  are  given  with  i|;e[0,TT).) 

THEOREM  3.15.  If  the  cylindroid  of  a  two-system  does 
not  lie  in  a  plane,  then  there  exists  a  unique  pair  of  lines 
on  the  cylindroid  that  intersect  at  right  angles. 

PROOF.  Let  <S-,>  and  <S2>  be  chosen  from  a  two-system 
such  that  their  central  axes  are  perpendicular  (that  there 
exists  a  pair  follows  from  the  Lemma).  Now,  all  pairs  of 
screws  with  perpendicular  central  axes  are  given  by  <c.S-,  + 
s^S2>  and  <-s.S^  +  c.S2>  where  t|;e[0,^).  By  Theorem  3.6, 
the  necessary  and  sufficient  condition  that  such  a  pair  of 
screws  also  have  perpendicular  central  axes  is 
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^""^^l    ■"  =+^2^  °  ^-^^^l    ^  ^^^2^  =  0 


Distributing  the  product  and  substituting  2h-,  =  S-.  o  S-i  , 

2h2  =  ! 

yields 


2h2  =  S2  o  S2,  d  =  -S^  o  S2 ,  33^  =  2s^c^,  and  c^^   =  cj  -  s^ 


(^2  -  ^^l)  ^2t  "  ^^2^  =  0 


Now,  it  is  not  possible  that  both  h2  -  h-,  =  0  and  d  =  0; 
for,  if  h-^  =  h2  and  d  =  0,  then,  by  Theorem  3.11,  all  the 
central  axes  would  lie  in  a  plane.  Thus,  there  is  a  unique 
solution  with  ^  e  [0,  — )  and  therefore  a  unique  pair  of 
screws  with  central  axes  that  intersect  at  right  angles. 
Q.E.D. 

The  two  screws  on  a  cylindroid  that  have  intersecting 
perpendicular  axes  are  called  the  principal  screws  of  a 
cylindroid.  This  pair,  of  course,  is  unique  for  cylindroids 
that  do  not  lie  in  a  plane.  If  the  cylindroid  does  lie  in  a 
plane,  then  any  pair  of  perpendicular  central  axes  can  be 
chosen  for  the  principal  screws. 

We  now  obtain  the  equation  of  the  cylindroid  and  the 
distribution  of  pitch  for  a  two-system.  Let  a  right-handed 
Cartesian  coordinate  system  be  located  such  that  the  x-  and 
y-axes  are  the  principal  axes  of  the  cylindroid  of  the  two- 
system  (thus,  the  z-axis  is  the  central  axis).  Let  <S-|>  and 
<S2>  be  the  principal  screws  where  s^  and  S2  determine  the 
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positive  X-  and  y-axis  directions.    Any  screw  <S>  in  the 

two-system  is  given  by  <S>  =  <c^S-]^  +  ^^^2^  where  ^p    e[0,  tt)  . 

The  position  of  the  central  axis  of  <S>  can  be  specified  by 

\|;,  which  is  the  angle  subtended  by  s  and  s-^,     and  by  the 

distance  z,  which  locates  the  intersection  of  the  central 

axis  with  the  z-axis  (see  Fig.  3.5). 

The  equation  of   the  cylindroid  can  be  obtained  by 

taking  the  reciprocal  product  of  S  with  the  rotor  S^  = 

-s.S-,  +  c  S^   whose  axis  is  perpendicular  to  both  the  z-axis 
V  1     i()  2 

and  the  central  axis  of  S  and  passes  through  the  origin. 
Using  eqn.  (3.6)  we  have 

S  o  s"^  =  (h  +  0)  cos  ^  -  z  sin  ^ 


=  z 


We  also  have 


S  o  S^  =  (c,S,  +  s,S-)  o  (-s.S?  +  cs^) 


=  -s ,c  (h.  -  h_) 
\\)    ^       I  2 


where  we  have  substituted  ^i  =  S  o  s, ,  h  =  S  o  S  ,  and  0  = 
S,  o  S-  =  S^  o  s  .  Thus,  the  equation  of  the  cylindroid  is 
given  by 
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Fig.  3.5.  Position  of  Central  Axis  on  Cvlindroid 
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z  =  -s^c^(h^  -  h2). 


The  distribution  of  pitch  is  given  by 


h  =  -J  S  o  S 


2^%^  ^  ^^2^  °  (%Sl  ^  %S2) 


2       2 

c,h-  +  s  ,h- 


where  we  have  substituted  h-,  =  —  S-i  o  Si  ,  h2  =  ^  So  o  Sp , 


and  0  =  S-j^  o  S2. 

Employing  the  identities   -s_  ,  =  s,c,,    -^(1  +  c^,)  = 

z   Z\\)         ^l   ^  2  24» 

2        1  2 

c  ,   and   —{1  -  c„,)  =  s,,   the  equation  of  the  cylindroid 
W  2  2^^  ■' 

and  the  distribution  of  pitch  can  be  expressed  as 


z  =  -  ^(h^  -  h2)s2^    , 


(3.7) 


h  =i(h^  -  h^)  ^|(h^  -  h2)C2^ 


(3.8) 


In  the  hz  plane,  it  is  evident  that  h  vs.  z  is  a 
circle.  Assuming  without  loss  of  generality  that  h-,  >  h2/  a 
circle  representation  of  the  two-system  is  illustrated  in 
Fig.  3.6.  For  h-j^  j^  h2/  each  point  on  the  circle  represents 
a  unique  screw  in  the  two-system;  for  h,  =  h2,  the  circle 
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Fig, 


6.  Circle  Representation  of  Two-System 
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degenerates  to  a  point,  which  represents  all  screws  in  the 
two-system  (the  cylindroid  lies  in  a  plane).  By 
consideration  of  the  circle  representation,  the  following 
facts  are  evident: 

(i)     The  principal  screws  assume  the  extreme  values  of 

pitch, 
(ii)    The  length  of  the  cylindroid  (i.e.,  along  the 

z-axis  )  is  |  h^^  -  h2  |  . 
(iii)   The  plane  containing  the  principal  axes  (i.e., 

the  X-  and  y-axes)  bisects  the  length. 


CHAPTER  IV 
ELEMENTARY  DYNAMICS 

In  this  chapter  we  express  basic  principles  of  dynamics 

in  terms  of  motors.   We  assume  familiarity  with  the  basic 

definitions  and  principles  of  analytical  mechanics.   We  only 

consider  finite  systems  of  mass  particles  (which  is  the 

realm   of   analytical   dynamics),   because   to   consider 

continuous  systems  would  introduce  unnecessary  complexities. 

Our  primary  concern  is  rigid  body  mechanics,   and  it   is 

sufficient  for  our  purposes  to  consider  that  a  rigid  body 

comprises  a  finite  number  of  particles. 

4 . 1   The  Wrench  and  Velocity  Motor  Revisited 
In  this  section  we  consider  some  additional  definitions 

and  properties  associated  with  the  wrench  and  the  velocity 

motor . 

DEFINITIONS.    A  pure   force   (rotation)   is  a  wrench 

(velocity  motor)  of  zero  pitch  and  a  couple  (translation)  is 

a  wrench  (velocity  motor)  of  infinite  pitch. 

We  have  that  any  wrench  can  be  expressed  as  a  sum  of 

pure  forces,  each  pure  force  of  which  corresponds  to  a  force 

in  the  force  system  of  the  wrench,  or  the  wrench  can  be 

decomposed  into  a  force  and  couple  (see  Section  3.3). 
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Before  we  can  give  physical  meaning  to  addition  of 
velocity  motors,  we  need  the  following 

THEOREM  4.1.  Let  V21  be  the  velocity  motor  of  frame  ^2 
relative  to  Z-,  and  V32  be  the  velocity  motor  of  frame  E^ 
relative  to  E2.   Then,  the  velocity  motor  of  £3  relative  to 

^1  is  V21  +  V32. 

PROOF.  This  follows  immediately  from  the  principle  of 
relative  velocities  for  points.   Q.E.D. 

Thus  any  velocity  motor  can  be  expressed  as  a  sum  of 
velocity  motors,  and,  in  particular,  it  can  be  decomposed 
into  a  rotation  and  translation  (see  Section  3.3). 

DEFINITION.  Two  force  systems  are  said  to  be 
equivalent  if  they  have  the  same  wrench.* 

All  the  common  reduction  procedures  for  equivalent 
force  systems  are  easily  deduced  using  motor  properties 
(e.g.,  that  forces  are  transmissible  and  that  any  two  forces 
acting  on  intersecting  lines  can  be  replaced  with  a  single 
force,  the  resultant,  on  a  line  through  the  point  of 
intersection) . 


*  This  is,  in  fact,  equivalent  to  the  standard  definition 
for  equivalent  force  systems  (the  term  equipollent  is  also 
sometimes  used).  In  dynamics  texts  the  definition  is 
usually  given  as  follows:  two  force  systems  are  said  to  be 
equivalent  if  the  sum  of  the  forces  for  each  system  are  the 
same  and  the  sum  of  the  moments  at  any  point  for  each  system 
are  the  same. 
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4  .  2   The  Momentum  Motor 

DEFINITION.  The  angular  momentum  distribution  of  a 
particle  system  is  called  a  momentum  motor.* 

Like  the  wrench,  the  momentum  motor  is  a  moment 
distribution  of  a  set  of  vectors,  namely,  of  the  linear 
momentum  vectors  (or,  simply,  the  linear  momenta)  of  the 
particle  system  (we  can  consider  these  momenta  point 
vectors).  Consider  a  system  of  mass  particles  with  masses 
m-i,  m2f  .  .  ,  mj^  located  at  points  ?-^,  P2  r  ..w  P^i  *  ^®*-  ^  ^^ 
directed  from  an  arbitrary  point  to  P^.  Then,  the  momentum 
motor  of  this  system  is  given  by 


V    i 
H  :   V   r   X  m.Vp 

i  =  l  i 


where  Vp   is  the  velocity  of  P^ .   The  principal  vector  is 
given  by 


K  =   y   m. v^ 
.  ^,   1  P. 
1  =  1      1 


which,  of  course,  is  the  linear  momentum  of  the  system. 
Denoting  the  center  of  mass  by  C,  the  velocity  of  the  center 
of  mass  by  Vp ,  and  the  total  mass  of  the  system  by  m  (thus. 


n 
m  =  ^      m . ) ,  we  have 
i=l   ^ 


*  It  appears  that  Clifford  (1878)  introduced  the  momentum 
motor;  he  called  it  the  momentum  of  twist.  Dimentberg 
(1965)  erroneously  states  that  Kotelnikov  (1895)  was  the 
first  to  define  the  momentum  motor. 
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h  =  mv^ 


We  say  that  the  momentum  motor  of  a  particle  system 
represents  the  total  momentum  of  the  system.  A  bivector  of 
the  momentum  motor  comprises  the  linear  momentum  and  angular 
momentum  (at  a  point)  of  the  system: 


(H)p  =  (mv^,  hp) 


4 . 3   Law  of  Momentum  for  a  Particle  System 
DEFINITION.   An  external  (internal)  wrench  acting  on  a 
particle  system  is  the  wrench  of  a  force  system  comprising 
of  external  (internal)  forces. 

We  shall  assume  the  strong  law  of  action  and  reaction; 
to  every  force,  there  is  an  equal  and  opposite  reaction 
force  with  the  same  line  of  action.*  As  is  customary  in 
elementary  dynamics  texts,  we  assume  that  this  is  equivalent 
to  Newton's  third  law.**  In  terms  of  motors,  we  have  that 
to  every  pure  force,  there  is  an  equal  and  opposite  reaction 


This   is   in  contrast  to  the  weak  law  of  action  and 


reaction,  which  states  only  that  forces  and  their  reaction 
forces  be  equal  and  opposite  (they  need  not  be  colinear). 

**  According  to  C.  Truesdell,  Newton  did  not  actually 
assert  this  law.  "In  Newton's  own  statement  of  his  third 
law,  there  is  no  explanation  of  what  kinds  of  'bodies'  he 
had  in  mind  or  what  he  meant  by  their  actions  on  each 
other,"  Essays  in  the  History  of  Mechanics,  Spr inger-Ver lag , 
New  York,  1968,  p.  27U. 
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pure  force.  Since  any  wrench  can  be  expressed  as  a  sum  of 
pure  forces  we  can  again  restate  this: 

NEWTON'S  THIRD  LAW.  To  every  wrench,  there  is  an  equal 
and  opposite  reaction  wrench. 

For  particle  systems  for  which  Newton's  third  law 
holds,  the  law  of  moment  of  momentum  also  holds.  If  H  is 
the  momentum  motor  of  a  system  and  M  the  external  wrench 
acting  on  the  system,  then  this  law  is  given  by 

M  =  H    ,  (4.1) 

where  the  derivative  is  taken  with  respect  to  an  inertial 
frame  (it  is,  of  course,  in  the  vector  field  form  in  which 
this  law  is  usually  given:  ra  =  -^rr)  .    Now,  if  f  denotes  the 

a  t 

net  force  of  the  force  system  and  mvp  the  total  linear 
momentum,  then  we  can  deduce  Newton's  second  law  by  simply 
equating  the  principal  vectors  of  eqn.  (4.1): 


f  =  mv^    .  (4.2) 


Newton's  second  law  is  also  known  as  the  law  of  linear 
momentum.  Since  the  law  of  linear  momentum  can  be  deduced 
from  the  law  of  moment  of  momentum  in  analytical  mechanics 
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(as  we  have  shown  here),*  we  call  eqn.  (4.1)  the  law  of 
momentum  for  a  particle  system. 

Equation  (4.1)  does  not  necessarily  characterize  the 
dynamics  of  a  particle  system;  that  is,  the  trajectories  of 
the  particles  cannot  in  general  be  determined  given  the 
external  wrench  as  a  function  of  time.  We  shall  not 
investigate  this  in  detail  since  it  is  with  rigid  bodies 
that  we  are  primarily  concerned,  and,  as  we  shall  show  in 
Chapter  IV,  eqn.  (4.1)  does  characterize  the  dynamics  for  a 
rigid  body.  Suffice  it  to  say,  that  for  an  arbitrary 
particle  system,  we  cannot  in  general  determine  the  particle 
velocities  at  an  instant  given  the  particle  positions  and 
the  momentum  motor  of  the  system  at  that  instant.  Consider, 
for  example,  any  two  particle  systems  for  which  the 
particles  have  equal  and  opposite  linear  momenta,  both 
vectors  of  which  are  parallel  to  the  line  on  which  the 
particles  lie.  The  momentum  motor  of  all  such  two  particle 
system  is  the  zero  motor. 

Of  course,  eqn.  (4.1)  can  be  applied  separately  to  each 
particle  in  a  particle  system,  in  which  case  the  motion  can 
be  completely  described.  In  fact,  we  need  only  apply 
Newton's   second   law   (i.e.,   eqn.  (4.2))   to  each  particle 


*    For  the  dynamics  of  deformable  bodies,  these  two  laws 
are,  in  general,  not  dependent. 
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(note  that  the  velocity  of  the  center  of  mass  of  a  single 
particle  system  is  simply  the  particle  velocity). 

4  .  4   Impulse,  Conservation  of  Momentum 
On  integrating  eqn .  (4.1)  over  a  time  interval  [t  ,  t] , 
we  obtain 

t 

/   MdT  =  H(t)  -  H(tQ)  =  AH  (4.3) 

t 
o 

t 

where   /   Mdx  is  the  impulse  motor  imparted  to  the  system. 

t 
o 

Thus,  we  have  that  the  impulse  motor  is  equal  to  the  change 
in  the  momentum  motor  AH.  If  M  =  0  then  H  will  remain 
constant  (i.e.,  AH  =  0);  this  is  the  law  of  conservation  of 
momentum  in  motor  form. 

We  say  that  the  impulse  motor  represents  the  total 
impulse  imparted  to  a  particle  system.  A  bivector  of  the 
impulse  motor  comprises  the  linear  impulse  and  the  angular 
impulse  imparted  to  the  system: 

t  t         t 

(  /   MdT)p  =  (  /   fdT,   /  mpdx) 

t^  t        t 

o  o        o 

where  f  is  the  net  force  acting  on  the  system.  With  ( H)  p  = 
(mv^,  hp) ,  it  follows  from  eqn.  (4.3)  that 


/   f d  T  =  mvp  I 

I-  ^   t 
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t  ^ 

/   mpd  T  =  hp  I ^     , 

t  o 

o 

from  which  we  can  deduce  the  conservation  laws  for  linear 
and  angular  momentum: 


mvc  1^   =  0 
o 


hp  It   =  ° 
o 


when  f  =  fflp  =  0. 

It  follows  that  the  total  momentum  is  conserved  if  and 
only  if  both  the  linear  momentum  and  angular  momentum  (at 
any  point;  hence,  all  points)  is  conserved. 

Finally,  we  note  that  if  a  constant  pure  force  acts  on 
a  particle  system,  then  the  angular  momentum  for  points  on 
the  central  axis  of  the  pure  force  is  conserved,  and  if  a 
couple  acts  on  the  system,  then  the  linear  momentum  is 
conserved. 

4 . 5  Power,  Work,  and  Energy  for  Rigid  Bodies 
Let  a  rigid  body  comprise  n  mass  particles  with  masses 
mj^  and  located  at  points  P^^ .  Let  V  be  the  velocity  motor  of 
the  rigid  body  frame  relative  to  some  inertial  frame  and  let 
M  be  the  wrench  of  the  force  system  acting  on  the  body.  vve 
assume  that  the  force  system  comprises  n  forces  fj^  (some  of 
which  may  be  zero)  acting  on  the  particles.   Then,  the  power 
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W,  which  is  the  time  rate  of  change  of  the  work  _W,  imparted 
to  the  body  by  the  force  system  is  given  by 

_W  =  It^     •  (Vy  +  V  X  OP^) 

=  (  I  f^)  •  Vq  +  V.  (  I  OP.  X  f.) 

=  M  o  V 

We  say  that  M  o  v  is  the  power  imparted  to  the  body  by  the 
wrench  M.  Now,  it  can  be  shown  (see  any  text  on  elementary 
dynamics)  that  for  rigid  bodies  the  internal  forces 
impart  no  power  to  the  body,  and  that  the  power  imparted 
to  the  body  (by  the  external  forces)  is  equal  to  the  time 

rate  of   change  of   kinetic  energy  I^,   which  is  given  by 

r  1      ? 
_T  =2,  ^  m^^Vp.  .   So,  we  have 


d   r  1    2 
^  =  dt  ^  2"i^P. 


Now,  let  H  be  the  momentum  motor  of  the  system.  Again 
choosing  an  arbitrary  fixed  point  0  in  the  body  frame,  we 
have  that 


I  ™iVpi  =  I(miVp.)  •  (Vy  +  V  X  OP^) 


=  (I  tniVp.  )  •  Vq  +  V  •  (  J^OPj^  X  (m^Vp.)) 
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=  K  •  Vy  +  V  •  h^ 


=  H  o  V 


Thus,  T  =  -^H  o  V  is  the  kinetic  energy  and 

M  o  V  =  ^(^H  o  V)    .  (4.4) 

Finally,  on  integrating  eqn .  (4.4)  over  a  time  interval 
[tQ,  t] ,  we  obtain  the  principle  of  work  and  energy: 


t 

t 


t  -,  ^ 

/       M    o    VdT    =       (-|h    o    V)     1^ 


t  t     ,  ^ 

where    /      MoVdx      (    =    /      WdT   =   W|)is      the     work      imparted 

t  to 

o  o 

over  the  time  interval. 

This  section  should  illustrate  the  central  importance 
of  the  reciprocal  product  as  regards  the  dynamics  of  rigid 
bodies,  as  this  product  is  intimately  related  to  the  concept 
of  power.  The  concept  of  power  is  important  tor  constrained 
systems:  if  the  total  power  imparted  to  a  system  by  the 
contraint  forces  is  always  zero  (i.e.,  tor  any  possible 
motion),  then  there  exists  a  formulation  of  the  equations  of 
motion  for  the  system  for  which  the  contraint  forces  are  not 
present.  Such  constraints  that  contribute  no  power 
(positive  or  negative)  to  the  system  are  known  as  workless 
constraints .  We  shall  give  the  definition  ot  workless 
constraints  in  terms  ot  motors  in  the  next  chapter. 


CHAPTER  V 
ELEMENTARY  KINETOSTATICS 

Phillips  (1984,  p.  4)  describes  the  term  kinetostat ics 

as  ".  .  .  the  study  of  angular  and  linear  velocities  and 

thus  the  kinematics  of  mechanism  on  the  one  hand,  this  being 

in  close  conjunction  with  the  study  of  forces  and  couples 

and   thus  of  massless  mechanism  on   the  other."    It   is 

important  to  point  out  that  it  is  instantaneous  kinematics 

to   which   Phillips   is   alluding   (which   is   the   study  of 

kinematic  properties  at  an  instant,  that  is,  the  relations 

of   velocities   and   higher  derivatives   in  mechanisms   and 

machines).    The   intimate   relation   between   instantaneous 

kinematics  and  statics  is,  perhaps,  best  illustrated  through 

the  use  of  motor  calculus.    We  consider  some  elementary 

aspects  of  these  subjects  and,  in  particular,  we  consider 

the   instantaneous   kinematics   and   statics   of   open-loop 

kinematic  chains  (or,  simply,  open-loop  chains). 

5.1   Constraint  Between  Bodies 


5.1.1   Characterizing  Constraint 

If  a  body  is  constrained  to  move  in  a  particular  manner 
relative  to  another  body,  then  we  shall  assume  that,  at  any 
instant,  the  constraint  between  the  bodies  is  completely 
characterized  by  a  motor  system  of  permissible  relative 
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velocity  motors  (Ball,  1900,  equivalently  characterized 
motion  freedom  via  screw  systems).  It  would  be  more  general 
to  consider  arbitrary  sets  of  velocity  motors,  as  opposed  to 
just  subspaces;  instantaneous  kinematic  properties  would 
then  become  much  more  involved.  This,  however,  is  beyond 
the  scope  of  this  work,  and  it  turns  out  that  in  many 
mechanisms  and  machines  (including  robots),  the  constraints 
can  be  sufficiently  modeled  by  motor  systems. 

We  shall  assume,  in  addition,  that  constraints  are 
workless,  which  we  now  define  in  terms  of  motors. 

DEFINITION.  The  constraint  between  a  body  1  and  body  2 
is  said  to  be  workless  if  M  o  V-|2  ~  0  where  M  is  the  wrench 
that  body  1  exerts  on  body  2  via  the  constraint  and  V-|2  ^s 
any  velocity  motor  of  body  1  relative  to  body  2  belonging  to 
the  motor  system  that  characterizes  the  constraint. 

We  note  that  -M  is  the  wrench  that  body  2  exerts  on 
body  1  (by  Newton's  third  law)  and  that  V2i  =  ~V-|2'  Since 
(-M)  o  V2T  =  Mo  V-,2'  ws  have  that  the  definition  is 
symmetric . 

We  call  the  wrench  exerted  via  a  workless  constraint  a 
constraint  wrench.    It  is  clear  that  a  constraint  wrench 


must  belong  to  a  reciprocal  motor  system.  It  is  thus 
possible  to  equivalently  characterize  a  constraint  by  a 
reciprocal  motor  system,  which  comprises  permissible 
constraint  wrenches. 
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5.1.2   Static  Equilibrium 

DEFINITION.  A  body  is  said  to  be  in  static  equilibrium 
if  the  body's  velocity  motor  with  respect  to  any  inertial 
frame  is  a  constant  translation  (or  the  zero  motor). 

It  follows  that  each  particle  in  the  body  has  a 
constant  velocity  so  that  the  net  force  on  each  particle  is 
zero.   Thus,  the  net  wrench  on  the  body  is  also  zero. 

DEFINITION.  A  body  is  said  to  be  grounded  (or,  it  is 
simply  called  ground )  if  it  is  in  static  equilibrium  and 
maintains  static  equilibrium  regardless  of  any  applied 
wrench.  A  body  is  said  to  be  constrained  to  ground  if  its 
motion  is  constrained  relative  to  a  grounded  body. 

Consider  a  rigid  body  constrained  to  ground  and 
stationary  relative  to  ground  (i.e.,  the  relative  velocity 
motor  is  zero).  Let  the  constraint  be  characterized  by  the 
motor  system  _N  and  let  M'  be  an  external  wrench  to  be 
applied  and  M^  be  the  resulting  constraint  wrench  also 
applied  to  the  body.  Since  the  body  is  initially  stationary 
relative  to  ground,  it  follows  that  the  body  is  in  static 
eqilibrium  before  M'  is  applied.  From  the  principle  of 
virtual  work,  we  have  that  the  necessary  and  sufficient 
condition  that  the  body  maintain  static  equilibrium  is 


M'  e  N^    ;  (5.1) 
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for,  this  is  equivilent  to  saying  that  M'  imparts  no  power 
during  a  so-called  virtual  displacement  (i.e.,  M'  o  v  =  0 
for  all  V  e  N) .  Since  the  net  wrench  exerted  on  a  body  must 
be  zero  if  the  body  is  in  static  equilibrium,  we  can 
determine  the  constraint  wrench  from 

M^  =  -M'    .  (5.2) 


5. 2   The  Screw  Pair 

We  shall  consider  the  screw  pair  illustrated  in 
Fig.  5.1  to  be  the  most  fundamental  joint  between  two 
bodies.  The  screw  pair  is  in  fact  the  most  general  single 
degree  of  freedom  lower  pair*  (Waldron,  1972).  Common 
multiple  degree  of  freedom  joints  such  as  the  cylindric, 
ball  and  socket,  planar,  and  universal  can  be  constructed  as 
special  serial  combinations  of  screw  pairs  (see  Duffy, 
1980)  . 

As  can  be  seen  in  the  figure,  uniform  rotation 
accompanies  uniform  translation:  a  rotation  of  9  about  the 
axis  of  the  screw  pair  accompanies  a  translation  of  he  along 
the  axis  where  h  is  the  pitch  of  the  screw  pair.  Should 
h  =  »,  however,  then  we  define  the  screw  pair  to  permit  only 


*    A  lower  pair  comprises  mating  surfaces  (of  the  joint 
between  two  bodies)  that  are  surfaces  of  revolution. 
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Fig.     5. 1.    The   Screw  P 


air 
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translation  for  which  the  pair  axis  can  be  taken  as  any  line 
parallel  to  the  translation  or  the  pair  axis  can  be  taken  as 
a  line  at  infinity.  For  h  ^  '^ ,  we  note  that  for  every  value 
of  9,  the  screw  pair  permits  a  screw  (i.e.,  a  one- 
dimensional  subspace)  of  velocity  motors  for  which  the  pitch 
is  also  h  and  the  central  axis  is  identical  to  the  axis  of 
the  screw  pair.  Since  the  screw  axis  is  fixed  in  the  two 
bodies  that  the  screw  pair  connects,  it  follows  that  the 
screw  is  fixed  in  both  bodies.  Thus,  the  screw  of  velocity 
motors  is  independent  of  the  value  of  9  (and  is  hence  the 
same  at  every  instant).  For  h  =  <»,  the  screw  of  velocity 
motors  at  any  instant  is  clearly  improper,  which  is  also 
fixed  in  both  bodies.  If  desired,  we  can  take  the  central 
axis  of  this  improper  screw,  which  is  a  line  at  infinity,  as 
the  central  axis  of  the  screw  pair. 

Two  practically  important  special  cases  of  the  screw 
pair  are  the  revolute  pair  for  which  h  =  0  and  the  prismatic 
pair  for  which  h  =  »,  both  joints  of  which  are  illustrated 
in  Fig .  5.2. 

Since  a  screw  pair  can  be  characterized  by  a  screw,  it 
can  be  characterized  by  any  motor  in  the  screw.  In 
particular,  it  is  very  convenient  to  represent  screw  pairs 
by  unit  motors.  Let  S  represent  a  screw  pair.  Then,  if  the 
pair  is  of  non-infinite  pitch  (as  in  Fig.  5.1),  then  we 
assume  the  direction  of  S  is  such  that  the  angle  of  rotation 
(9  in  Fiy.  5.1)  is  measured  in  a  right  sense  about  s;  it  the 
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5  1 


(a) 


^ 


(b) 


Fig.  5.2.  The  (a)  Revolute  and  (b)  Prismatic  Pair 
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pair  is  prismatic  (as  in  Fig.  5.2b),  then  we  assume  the 
direction  of  S  is  such  that  the  length  of  translation  ( I  in 
Fig.  5.2b)  is  measured  in  a  positive  sense  in  the  direction 
of  3.  The  velocity  motor  can  then  be  given  by  qS  where  q  is 
the  measure  of  displacement  that  is,  the  joint  displacement 
of  the  screw  pair  (q  =  9  in  Fig.  5.1  and  q  =  I  in 
Fig.  5.2b). 

Finally,  we  note  that  by  Theorem  3.5,  a  unit  motor  that 
represents  a  screw  pair  is  fixed  in  both  bodies  connected  by 
the  pair. 

5 . 3   Open-Loop  Chains 

5.3.1  Modeling 

We  shall  assume  that  the  joints  connecting  successive 
bodies  (or  links)  in  an  open-loop  chain  are  screw  pairs. 

Let  an  open-loop  chain  comprise  n  links  numbered  from 
the  base  outward  as  illustrated  in  Fig.  5.3.  We  number  the 
base  (or  the  ground  link)  0.  The  joints  are  denoted  by  the 
unit  motors  S^  where  the  i  joint  connects  links  i-1  and  i, 
and  the  corresponding  joint  displacement  is  given  by  q^^ 
where  q^^  =  B  ^    {q^    =  Z^)     if  Sj^  is  proper  (improper). 

5.3.2  Velocity  and  Acceleration  Analysis 

Obtaining  the  equations  that  are  necessary  for  velocity 
and  acceleration  analysis  of  open-loop  chains  is 
straightforward  using  motor  calculus. 
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Fig.    5.3.    The   Open-Loop    Chain 
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Let  V4  denote  the  velocity  motor  of  the  j^^  link 
relative  to  the  ground  link  (thus,  Vq  =  0).  Then,  by 
Theorem  4.1,  we  have 


V^  =   t   S.q.    , 
3         ■  ^,   11 


(5.2) 
i  =  l 


which  can  be  expressed  in  the  recursive  form 

Vj  =  Vj_;L  +  Sjq^    .  (5.3) 

Differentiating  both  sides  of  eqn.  (5.2)  (with  respect 
to  the  ground  frame),  we  obtain 


V.  =   y   (S.  q.  ) 

J    i=l    ^  ^ 


Now,  since  Sj^  is  fixed  in  both  links  i-1  and  i,  we  have  that 

S-  =  V._,  X  s.  =  V.  X  s.   (the  time  derivatives  of  Sj^  with 

respect   to   links   i-1   and   i   are  both  zero)  ;   we  choose 

S.  =  V.  X  s.   for  the  following.   Thus,  we  have 
111  ^  ' 


i 

V.  =   y   (S.  q.  +  V.  X  S.  q. 
J     . ^,     1^1     1     1^1 


=   )   [S.q.  +(y   S,  q,  )><S.q. 
i  =  l    ^     '^  k=l  ■'■ 
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=   fs.q.+   y    TS,   xq.q    .  5.4 

1=1         1=1  k=l 


Differentiating  both  sides  of  eqn.   (5.3)   we  obtain  the 
recursive  form 


V.  =  V.  ,  +  S.q.  +  V.  X  S.q     *  (5.5) 


SuiTunarizing  , 


h  =  A  ^i^i   ' 


V^    =       I      S,q;  ,  (5.2) 

i  =  l 


V.=       I       S.q.+       y  y      S,      xS.q.q,  ,  5.4 

J         .  '"t       1^1         .  ^T    ,  ^T       k  i^i^k 

-•         1  =  1  1  =  1    k  =  l 

which  are  explicit  expressions  in  the  first  and  second  time 
derivatives  of  the  joint  displacements,  and 


Vj  =  W^_^    +  S^q.    ,  (5.3) 


V.  =V,+Sq.+V.xS.q.  (5.5) 


which  are  recursive  expressions. 
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5.3.3   Statics 

We  shall  assume  that  open-loop  chains  are  fully 
actuated;  that  is,  associated  with  each  joint  is  an 
actuator.  (Most  industrial  robots  would  fall  into  this 
category.)  Each  actuator  in  a  fully  actuated  open-chain  can 
create  a  wrench  "in-between"  the  links  joined  by  the 
associated  joint.  We  now  make  these  notions  precise  and 
obtain  the  equations  of  equilibrium  for  the  open-loop  chain. 

We  denote  the  screw  of  the  i  actuator  by  S.  and  we 
assume  that  this  screw  is  fixed  in  the  i  link.  The  wrench 
that  the  i  actuator  exerts  on  the  i  link  is  given  by 
^i^i'  ^"^  *-^®  wrench  that  is  exerted  on  the  (i-1)  link  is 
~^i^i  ('^s'^ce,  the  actuator  creates  a  wrench  "in-between"  the 
two  links).  We  note  that  x^  is  the  magnitude  of  a  force  or 
torque  (which  we  allow  to  be  negative)  depending  on  whether 
S^  is  proper  or  improper. 

Suppose  that  on  each  link  of  an  n-link  open-loop  chain 
an  external  wrench  M[  is  exerted  and  that  we  require  the 
chain  be  in  static  equilibrium;  that  is,  each  link  be  in 
static  equilibrium  (we  assume,  of  course,  that  the  base  is 
grounded )  . 

Then,  applying  condition  (5.1)  to  link  i,  we  obtain 


T  ■  s?  + 


n 

y  M   e  <S.  >' 

^  .  p      1 

p=i  ^ 
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and  thus  the  constraint  wrench  M^  e  <S^>   is  determined  from 


n 
M^  =  -i\Sl    +       I      H')  .  (5.6) 

p=i   ^ 

n 
Now,  Sj^  /i     <S^>^,     for,  if  sf  £  <S^>'^  and   I      n'    j6    <S.>  , 

n  p=i   ^ 

then  1      M'    could  not  be  equilibrated  (i.e.,  eqn.   (5.6) 

p=i   P 
would  be  inconsistent).    Taking  the  reciprocal  product  of 

both  sides  of  eqn.  (5.6)  with  S  ^^ ,  we  obtain 


n 

:.Sf  o  S.  +  (  y   M'  )  o  S.  =  M^  o  S. 
D  J     1     p=i   P      1     1     ^ 


Since  S^  o  Sj^  /^  0  and  M^  o  s^  =  0 ,  we  obtain 


n 

(  I.  M')oS. 

T  =  -  P=\      .  (5.7: 

1   1 


We  thus  have  a  relation  between  the  actuator  wrench 
magnitudes  t-,,  12/  .../  t^^  and  the  applied  wrenches  M]  ,  M^/ 
...,  Mjlj .   If  desired,  the  constraint  wrenches  M^ ,  M2 ,  ..., 


M^  can  be  subsequently  determined  from  eqn.  (5.6). 

In  industrial  robots,  it  is  very  coimnon  that  joints  be 
either  revolute  or  prismatic  pairs.  If  S:  represents  a 
revolute  pair  (thus,  S^^  is  a  rotor),  then  we  shall  always 
assume 


96 


sf  =  sr 


if  Sj^  represents  a  prismatic  pair  (thus,  S^  is  improper), 
then  we  shall  always  assume  S^  is  a  rotor  such  that 

(Sf)"  =  Si    . 

It  is  in  this  manner  that  industrial  robots  are  usually 
actuated.  We  note  that  if  S^  represents  a  revolute  or 
prismatic  pair,  then  5^03^=  1,  and  eqn.  (5.7)  simplifies 
to 

n 

J.    =    -{     I      M|)  o  S.    .  (5.8) 

i  =  j         ■> 

Equation  (5.7)  is  the  important  result  of  this 
subsection.  By  using  D'Alembert's  principle,  we  can  extend 
this  equation  so  that  a  complete  dynamic  analysis  can  be 
performed  for  an  open-loop  chain.  We  consider  this  in  the 
next  chapter. 


CHAPTER  VI 
DYNAMICS  OF  RIGID  BODIES 

In  this  chapter  we  tie  together  many  of  the  results  of 
previous  chapters  in  order  to  derive  the  equations  of  motion 
for  rigid  bodies  in  terms  of  motors. 

We  first  give  a  novel  derivation  of  the  inertia 
dyadic  (the  inertia  dyadic  was  introduced  by  von  Mises, 
1924a).  Moreover,  by  using  a  dyad  expansion,  we  obtain  a 
particularly  simple  form  for  this  dyadic. 

We  then  obtain  the  equations  of  motion  for  a  rigid 
body;  these  equations  are  simply  specializations  of  the 
equations  for  a  particle  system,  which  were  given  in 
Chapter  IV. 

Finally,  the  equations  of  motion  for  an  actuated  open- 
loop  chain  are  derived.  It  is  believed  that  these  equations 
are  simpler  in  form  than  any  others  found  in  the  literature. 

6. 1   The  Inertia  Dyadic 
6.1.1   Derivation 

Let  us  consider  a  rigid  body  of  mass  m  and  let  I^  be 

the   inertia   tensor   about   the   mass   center   C.  (The 

superscript   in   l^   indicates   that   it   is   a  linear 

transformation   for   vectors.)     In   addition,   let  r^      be 

directed  from  an  arbitrary  point  to  the  mass  center.  If  the 
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body  has  velocity  motor  V  (with  respect  to  an  intertial 
frame),  then  the  momentum  motor  is  given  by 


H  :  r^,  X  mV(.  +  I^v  (6.1) 


where  K  =  mvp. 

Let  p-^,  P2,  P3  be  unit  vectors  each  of  which  is  in  the 
direction  of  a  principal  axis  (thus,  these  are  the 
eigenvectors  of  I^).  With  r-j^,  r2,  r^  denoting  the 
respective  radii  of  gyration  we  have  the  following  spectral 
decomposition : 


I^  =  m   I   r^p.p. 


where  we  assume  that  standard  definition  for  vector  dyads 

(i.e.,  (ab)c  =  a(b  •  c) ) .   We  substitute  this  dyadic  and  the 

3 
identity  dyadic    I      p.p.  into  eqn.  (6.1),  which  yields 

i=l   ^  ^ 

3  ^    2      ~ 

H  :  rp  X  m(  I      p.p.  )v   +  m(  j;   r  p  p  )v 

i=l        ^      i=l   ^  ^  ^ 


(The  motivation  for  introducing  these  vector  dyadics  will  be 
soon  apparent.)  Exploiting  properties  of  vector  dyads,  we 
obtain 


?  ?    2  ^ 

H  :  m  I      (r   X  p. ) ( p.  •  v  )  +  m  I      r  p  (p   •  v) 

i=l  ^       i=l 


(6.2) 
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Now,    we    define    P^    and    P°°   such    that 


Pi     :    r^    X    p.  , 


P.     :    p.  , 

1         *^i 


and  we  note  that 


v^  =  P.  o  V 

C     1 


V  =  P.  o  V 

1 


Substituting  into  eqn.  (6.2),  we  obtain 


3  3 

H  =  m   y   PrCP'r  o  V)  +  m   y   rl"  P.  (P.  o  v) 
1=1  1=1 


.  (6.3) 


Thus,  H  is  a  linear  function  of  V.   Defining  the  motor 
dyadics , 


lO  =  m 


0^0 


.  '- ,   1  1 


1  =  1 
3 


2_,"„'» 


I-  =  m  I   r1^p"p: 

i=l   '  '  ' 


(6.4) 


(6.5) 


r*        no        9   00   oo 

I  =  m   y    P*:^P^  +  r^^P.  P.  ) 
i=l    ^    '  '    '    ' 


(6.6) 


we  have  that 
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H  =  IV  (6.7) 


where 


I  =  1°  +  I"    . 


We  call  I  the  inertia  dyadic  of  the  body.   We  note  that 


IV  :  r^  X  mV(^    ,  (6.8 


I^V  :  I^v    .  (6.9 


It  is  clear  that  I^  is  determined  solely  by  the  body's  mass 
and  the  mass  center  and  that  l"  is  determined  solely  by  the 
mass  distribution. 

Finally,  we  note  that  {p{ ,  pO ,  p^ ,  p^,  p^,  pp  is  a 
standard  basis  fixed  in  the  body  frame.  It  is  clear  that 
the  principal  axes  of  inertia  constitute  the  Cartesian 
coordinate  system  associated  with  this  basis  (i.e.,  the 
central  axes  of  P^ ,  p9 ,  p9 )  . 

6.1.2   Dyad  Expansion 

A  dyad  expansion  for  an  inertia  dyadic  I  can  be  rapidly 
determined  by  using  the  following  identity: 
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3 

iil    '    '  '    '    ^ 


,    0  00        0  00 

(P.  -  r.P.)(P^  -  r.P.)]    .  (6.10) 


Introducing  the  motors 

Pi  =  P?  ^  r,p-    ,    P2  =  P?  -  -iPr    ' 

P3  =  pO  +  ^2^2    '    P4  =  P2  -  ^^2^2    ' 

P5  =  P§  -  r3Pj    ,    P6  =  P§  -  -3^1    ' 
into  identity  (6.10)  we  can  then  rewrite  definition  (6.6)  as 


6 

I  =  f   I   P^P.    .  (6.11) 

^  i=l   ^  ^ 


It   is   easily   verified   that   the   set   {P- }   is   co- 
reciprocal  and  that 


Pj  o  P^  =  2hj 


where  h-  is  the  pitch  of  P^  and 
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h^  =  r^    ,    h2  =  -r^ 


h3  =  r2    ,         h^  -  -r 2 


h^  =  rj    ,    hg r3 


In  addition,  {Pj^}  consists  of  eigenmotors ,  for, 


6 

,m 


1  =  1        -* 


=  fp.(P,  o  P^) 


=  mh. p. 
J  J 


the  corresponding  eigenvalues  are 


Xt  =  mr-j^  (-  mh-]^)    ,    X2  =  -mr -^    (=  mh2)    r 


X-j  =  mr2  {=  mh^)    /    X^  =  -mr2  (=  mh^)    , 


Xr  =  mr-^  (=  mhc)    ,    Xg  =  -mr^  (=  ™hg) 
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We  have  tacitly  assumed  and  we  shall  continue  to  assume 
that  the  radii  of  gyration  (i.e.,  r-,,  r2,  r-,)  are  all  non- 
zero.  It  follows  from  this  assumption  that  (P^  }  is  a  basis 

of  eigenmotors  and  that  eqn.  (6.11)  gives  the  desired  dyad 

X. 

expansion  in  which  the  simplification  :=: — =—  =  t:    has  been 

p.  op.    2 
1   1 

incorporated . * 

The  screws  <Pi>  ai^e  what  Ball  called  the  principal 
screws  of  inertia  of  a  rigid  body  (Ball  actually  used  the 
synthetic  definition  for  the  screw).  It  is  clear  that  the 
motors  Pj^  are  fixed  in  the  body. 

Finally,  using  the  dyad  expansion  for  I,  we  can 
immediately  write  down  the  dyad  expansion  for  its  inverse 
(see  subsection  2.5.3): 

,     6    X-1 

^     =  J,  p-k-  ^^i 

1  =  1   1   1 


1    ^   1 

^   I   ^  P.  P.    .  (6.12) 

2m  .^^  j^2   1  1 


*  Equation  (6.11)  would  still  be  correct  even  if  some  of 
the  r^  were  zero  because  identity  (6.10)  would  still  hold. 
It  would  not,  however,  be  a  dyad  expansion  according  to  our 
definition.   For  example,  suppose  r^^  =  0  and  ^2    /^  0 ,  r^  ^ 

0.   Then,  P^  =  P2  (=  P^)  and  eqn.  (6.11)  can  be  rewritten  as 

6 
I  =  mP^P^  ^  H   5;   p.p.    . 

1  =  3 
The  set  (Pj^}  still  comprises  eigenmotors,  but  I  cannot  be 
expressed  as  a  dyad  expansion  since  P-.  0  P-i  =  0. 


104 


6.1.3   The  Intertia  Inner  Product 

We  now  obtain  a  particularly  useful  inner  product  as 
will  be  demonstrated  in  section  6.3  where  we  obtain  the 
equations  of  motion  for  an  open-loop  chain.  This  inner 
product  is  induced  by  an  inertia  dyadic,  which  we  call  the 
inertia  inner  product. 

We  show  that  an  inertia  dyadic  I  does,  in  fact,  induce 
an  inner  product  by  using  the  dyad  expansion: 


m         6 
I   =  ?    I    P.P. 

2  iii    ^  ^ 


Now,  we  must  show  that  I  is  symmetric  and  that  M  o  IM  >  0  if 

and  only  if  M  7^  0.   Since  each  dyad  P^P^  is  symmetric,  it  is 

clear  that  I  is  symmetric.    In  order  to  demonstrate  the 

positive  def initeness ,  let  some  motor  M  be  expressed  as  a 

6 
linear  combination  of  the  motors  P^  ,  that  is,  M  =   T   a. P.  , 

1=1 
and  then  we  have 


m    6  -> 

MolM  =  ^  I       (P.  oM) 

^  i  =  l    ^ 


!  j^  '^  °  "i''  "i 


6      ~    ^ 
y   (mhf^)  a'^ 
1  =  1 
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Since  all  h-  j^    0,  we  have  that  \      (mh.)  a.  >  0   if  and  only 

i  =  l 
if  not  all  Oj^  =  0.   Since  not  all  a^^  =  0  is  equivalent  to 

M  =  0,  we  have  that  M  o  IM  >  0  if  and  only  if  M  ;^  0;  thus,  I 

induces  an  inner  product  (  • |  •)  . 

A  useful  expression  for  (M|N)  can  be  obtained  using 

eqns.  (6.8)  and  (6.9).   With  C  denoting  the  mass  center  of 

the  body  to  which  I  corresponds,  we  have 


M   o    I^N  =    (M)^    o    (  I°N)(;> 


=    ( m ,    nip  )     o    { min ,    0 ) 


=   mnip    o   iip 


and 


M   o    i"n   =   m    •    I'^n 


thus. 


(MlN)    =    M    o    (lO    +    I~)N 


=    M    o    l^N   +    M    o    i"n 


=   mnip    •    lip   +   m    •    In 


(6.13) 
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We  call  the  norm  induced  by  an  inertia  inner  product 
and  inertia  norm.*  We  can  give  an  interesting  physical 
interpretation  for  the  inertia  norm.  From  Chapter  IV,  we 
have  that  the  kinetic  energy  is  given  by 


jr  =  ^  H  o  V 


Now,  since  H  =  IV  we  have 


I^  =  ^V  o  IV 


=  ^  «VII^    .  (6.14) 


Finally,  using,  eqn.  (6.13),  we  obtain  the  common  expression 


X=  ^v^  *  '^C  ^  i^  *  ""-^^    *  (6.15) 


6 . 2   Equations  of  Motion  for  a  Rigid  Body 
Before  we  derive  the  equations  of  motion  for  a  rigid 

body,  we  need  the  following: 

DEFINITION.    A  dyadic  A  is  fixed  in  a  frame  if  for 

every  fixed  motor  M,  the  motor  AM  is  also  fixed. 


*  In  some  advanced  dynamic  texts,  the  metric  d( •, •)  such 
that  d(M,N)  =  DM-Nll,  where  II  •  II  is  an  inertia  norm,  is  called 
the  inertia  metric. 
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It  is  clear  that  if  the  dyads  that  constitute  a  dyadic 
comprise  fixed  motors,  then  the  dyadic  is  fixed.  In 
particular,  we  have  that  the  inertia  dyadic  of  a  body  is 
fixed  in  the  body  frame. 

Consider  a  rigid  body  for  which  I  is  the  inertia 
dyadic,  V  is  the  velocity  motor,  and  M  the  applied  external 
wrench.  Then  substituting  H  =  IV  into  the  law  of  momentum 
(eqn.  (4.1)),  we  obtain 


"'ft" 


I^IV  +  V  X  IV  (6.16) 


d' 
where  -r—    denotes  differentiation  with  respect  to  the  body 

frame.   Since  I  is  fixed  in  the  body  frame,  we  have 


d'       d' 
IV  =  I V 

dt^^    ^dt^ 


But 


d' 

-—V  =  V  -  V  X  V 

dt 


=  V 


d ' 
so  that   tttIV  =  rv,   and  we  can  write  the  law  of  momentum 

for  a  rigid  body  as 
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M  =  IV  +  V  X  IV 


(6.17) 


The  law  of  conservation  of  momentum  for  rigid  bodies  is 
given  by 


(IV) 


=  0 


(6.18) 


where  M  =  0  over  the  interval  [t^,  t]  ,  and  the  principle  of 
work  and  energy  for  rigid  bodies  is  given  by 


/    M  o  VdT  = 


(^nvi^) 


(6.19) 


We  note  that  by  using  motor  calculus,  the  equations  of 
motion  for  a  rigid  body  attain  a  very  simple  form.  The 
above  dynamical  laws  are  essentially  in  the  same  form  as 
given  by  von  Mises  (1924a  and  1924b)  (although,  von  Mises 
did  not  employ  an  inertia  norm);  our  derivations,  however, 
are  quite  different. 


6 . 3  Dynamics  of  an  Open-Loop  Chain 
The  equations  of  motion  for  an  open-loop  chain  are 
readily  derived  using  D'Alembert's  Principle  and 
eqn.  (5.7).  We  use  the  same  notation  as  in  Chapter  V,  and, 
in  addition,  we  denote  the  momentum  motor,  inertia  tensor, 
and  energy  inner  product  tor  the  i  link  by  H^ ,  I^^,  and 
(•!•)..   The  intertial  wrench  "applied"  to  the  i*-^  link  is 
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given  by  -H.  where  H^  =  I^Vj^.  We  obtain  the  requisite 
actuator  torque/force  magnitudes  by  extending  eqn.  (5.7) 
(i.e.,  by  using  D'Alerabert's  Principle): 


n 

•i  =  -^i[    I      (Mp  -  Hp)]  °  S.  (6.20) 

p=i 


where  k^  =  (S^  o  Sj^  )"■'•.   The  constraint  wrenches  are  then 
computed  from  from 

M^  =  H.  -  (M|  +  T.S.  )    .  (6.21) 


Since 


n  n 

•■  =  -k.  (  y   M'  )  o  S.  +  k.  (  y   H  )  o  S. 

1      i^.p  1     i^.p      1 

p=i  ^  p=i  ^ 


it  is  clear  that  we  can  compute  components  of  actuator 

torque/force  magnitudes  due  to  the  MV  and  due  to  the  -H 

P  p 

separately.    We  thus  restrict  ourselves  to  the  dynamic 

problem;  that  is,  we  assume  M,l  =  0  for  i  =  l,  2,  ...,n.   In 

P 

addition,  we  define 


T  . 

^i  =  ^  (6.22) 

i 


(note  that  t[  =  Tj^  for  open-loop  chains  whose  joints 
comprise  only  revolute  and  prismatic  pairs).  Then, 
eqn.  (6.20)  reduces  to 
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n 

:\    =    {     I      H    )    o    S.        .  (6.23) 

P=i 


We  now  "expand"  the  right  hand  side  of  eqn.  (6.23)  in 
order  to  obtain  the  so-called  explicit  form,  which  is  a  form 
explicit  in  the  q.  and  q.  .  Substituting  for  H  in 
eqn.  (6.23)  and  distributing  the  reciprocal  product,  we 
obtain,  ■ 


T!  =  [y   (IV   +V   xiv)]oS. 
1     >.PP     P     PP       1 


p=i 


n 

=   7(3.  oIV+S.  oVxIV) 
p=i    1     P  P     1     P     P  P 

Employing  the  identity  S^  x  v^  o  T^V,,  =  S.:  o  v^  x  t  v^  anc 

i     P     PP     ■•■     p     PP 

introducing  the  inertia  inner  product  (  •!  •)    yields 

n 

t!  =   y   [(S.  I V  )   +  (S.  X  V  I V  )  ] 
1    p^i     1  '  P  P      1     P'  P  P 


Substituting  for  V   and  V    using  eqns .  (5.2)  and  (5.4)  and 
rearranging,  we  obtain 


n         p     ..      p    J 
T-  =  I       [(S.  I  I       S  q   +   1    I   S   X  S  q  4,,)„ 
^    p=i     ^  j=l   J  ^    j  =  l  k=l   "^     J  J  K  P 


(S.  X   I       S.q.  I  y   S,  q,  ) 
j  =  l   -•  -'  k=l 


Ill 


pi  'X  "^'^^'^'^ 


I    i   (S.  |S   X  S  )q  qk 
j=l  k=l    ^   '^     J   ^ 


'  jl  Jl  '^'    ^    ^jl^k^'^j^k^     • 


(6.24) 


The  equations  of  motion  can  thus  be  expressed  in  the  form 


n    n 


^1  =  Jl  '^i'^  *  ,k  Jl  '^^^'^^^     • 


(6.25) 


The  coefficients  D^^  and  C^.|^,  determined  from  eqn .  (6.24), 
can  be  expressed  in  the  form 


n 


ij 


p=max( i , j ) 


(6.26) 


and 


n 


Cllk  =   I   KS^  "    S.|S,)p  -  t^,(Si|S,  «  S,)pl   (6.27) 
p=max( i, j ,k) 


where 


-jk 


We  can  obtain  an  even  simpler  form  for  these  coefficients  by 

* 
introducing  an  inner  product  ('I'^q   defined  as  follows: 
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(M|N)   =  I       (M|N)      .  (6.28) 

p=q 

This  new  inner  product  is  actually  an  inertia  inner  product 

n 
induced    by    the    composite    inertia   dyadic     1      I  . 

p=q   P 
Substituting  definition  (6.28)  into  the  expressions  (6.26) 

and  (6.27),  we  obtain 

D..  =  (S^|S. )     ,    p=max(i,j)  (6.29) 

and 

C..,  =(S.  xSlS,  )   -t.,(S.  Is.  xS.  )*   ,     (6.30) 
ijk      1     j'  k  p     jk   i'  J     k  p 

p=max( i , J ,k) 


Featherstone  (1984,  pp.  80-81)  obtained  essentially  the 
same  form  for  D^^  given  by  eqn.  (6.29)  in  his  "Composite- 
Rigid-Body  Method."  It  is  believed,  however,  that  the  form 
for  Cj^-jj^  given  by  eqn.  (6.30)  is  new. 

Finally,  we  note  that  since  C-  ^^  is  the  j,k  entry  of 
the  matrix  of  a  quadratic  form,  we  can  replace  it  with  C|^|^ 
where  for  j  >  k, 


c:  ,  =  -^(C.  .,  +  C,  .  ) 
ijk    2   ijk     ikj 


1  *  * 

=  -^[(S.  X  Sis.)   +  (S,   X  s  |S.) 
2i     jkp      k     j'lp 


(S.  X  s,  |S  )  1    ,    p  =  max(i,j,k)    , 
1     1^   J  P 
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and  for  j  <  k. 


C    =  C 
ijk     ikj 


6.4   Discussion 


The  notion  of  coordinate  independence  is  of  central 
importance  in  this  chapter.  We  did  not  introduce  a 
coordinate  system  (of  motors)  in  order  to  construct  the 
inertia  dyadic,  nor  did  we  introduce  any  coordinate  systems 
in  order  to  derive  the  equations  of  motion  for  an  open-loop 
chain.  It  is  believed  that  these  coordinate  independent 
derivations  are  new. 

Algorithms  for  dynamics  computations  can  be  developed 
once  coordinate  systems  are  introduced.  Ball  (1900) 
expressed  dynamic  equations  of  a  rigid  body  in  terms  of 
coordinates  referred  to  a  basis  of  eigenmotors  of  the  body's 
inertia  dyadic  (the  P^  in  definition  (6.1)).  Most  other 
researchers  who  have  applied  motor  calculus  or  screw  theory 
to  rigid  body  dynamics  have  used  standard  bases  (e.g., 
von  Mises,  1924a  and  1924b;  Dimentberg,  1965;  Yang,  1971; 
Woo  and  Freudenstein,  1971;  and  Featherstone ,  1983  and 
1984)  . 

In  order  to  develop  efficient  dynamics  algorithms  for 
open-loop  chains,  great  care  must  be  taken  in  the  selection 
of   both   coordinate   systems   and   the   manner   in   which 
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coordinates  are  transformed  from  one  system  to  another, 
Featherstone  (1984)  has  investigated  this  subject 
extensively. 

There  has  been  a  vast  amount  of  literature  recently  as 
regards  the  development  of  efficient  dynamics  algorithms  for 
open-loop  chains  using  vector  methods  (see,  for  example, 
Orin,  McGhee,  Vukobratovic ,  and  Hartoch,  1979;  Luh,  Walker, 
and  Paul,  1980;  Walker  and  Orin,  1982;  Thomas  and  Tesar, 
1982;  Renaud,  1984;  Wang  and  Kohli,  1985;  and  Lathrop, 
1985).  Only  Featherstone  (1984)  has  considered  the  use  of 
motor  calculus  in  developing  efficient  algorithms. 

A  comparison  of  the  forms  of  the  equations  of  motion 
for  an  open-loop  chain  given  in  this  chapter  and  in 
(Featherstone,  1984)  with  the  equations  given  in  the  above 
cited  references  that  use  vector  methods  would  illustrate 
the  greatly  simplified  structure  that  results  by  using  motor 
calculus.  Perhaps  an  important  contribution  of  this  chapter 
is  the  simple  forms  obtained  for  the  C^-j^  coefficients.* 
The  simplicity  can  be  realized  by  comparing  with  the  forms 
obtained  in  (Thomas  and  Tesar,  1982)  and  (Wang  and  Kohli, 
1985).  In  both  these  papers,  elaborate  tabular  schemes  are 
necessary  in  order  to  determine  these  coefficients.    The 


*  Featherstone  (1984)  did  not  determine  expressions  tor 
these  coefficients,  as  the  algorithms  that  he  developed  did 
not  require  dynamic  equations  explicit  in  the  joint 
velocities. 
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forms  given  in  this  chapter  require  no  such  schemes;  the 
tables  are,  in  effect,  incorporated  in  the  algebra. 


CHAPTER  VII 
CONCLUSIONS 

Motor  calculus  is  clearly  well-suited  for  rigid  body 

dynamics.   This  has  been  demonstrated  by  von  Mises  (1924a 

and  1924b)  and  by  several  recent  contributions  (perhaps, 

most  notably,  Dimentberg ,  1965  and  Featherstone ,  1983  and 

1984).    A  main  objective  of  this  dissertation  has  been  to 

treat  motor  calculus  and   its  application  to  rigid  body 

dynamics  in  a  more  general  and  systematic  manner. 

Important  contributions  include: 

the  new  motor  definition, 

the  development  of  a  dyad  expansion  for  a  motor 

dyadic , 

a  rigorous  treatment  of  motor  derivatives  and 

integrals , 

a  simple  derivation  of  the  circle  representation 

for  a  two-system, 

simple  derivations   and   representations  ot   the 

inertia  dyadic  of  a  rigid  body  and  the  equations 

ot  motion  for  an  open-loop  chain. 

It   is   hoped   that   some   ot   the   foundational   areas 

presented   in   this  dissertation  will   facilitate  the  many 

potential  applications  ot  motor  calculus.   Two  suggestions 

for  future  work  are  the  investigation  ot  the  dynamics  of 
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general   mechanisms   and   further   research   on   developing 
efficient  algorithms  for  dynamic  analyses. 
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